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ABSTRACT 

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation 
(6.0eU < hu < 13.6eU) from hot early type OB stars on clumps in star- forming molecular 
clouds. Clumps in FUV-illuminated regions (or photodissociation regions or PDRs) undergo 
external heating and photodissociation as they are exposed to the FUV field, resulting in a 
loss of cold, molecular clump mass as it is converted to warm atomic gas. The heating, if 
rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks 
into the clumps, compressing them to high densities. The clumps lose mass on relatively 
short timescales. The evolution of an individual clump is found to be sensitive to three 
dimensionless parameters: the ratio of the initial column density of the clump to the 

column N 0 ~ 10 21 cm~ 1 2 of a warm FUV-heated surface region; u, the ratio of the sound speed 
in the heated surface to that, in the cold clump material; and t FUV /t c , the ratio of the “turn-on 
time” t FUV of the heating flux on a clump to its initial sound crossing-time t e . The evolution 
also depends on whether a confining interclump medium exists, or whether the interclump 
region has negligible pressure, as is the case for turbulence-generated clumps. In this paper, 
we use spherical 1-D numerical hydrodynamic models as well as approximate analytical 
models to study the dependence of clump photoevaporation on the physical parameters of 
the clump, and to derive the dynamical evolution, mass loss rates and photoevaporative 
timescales of a clump for a variety of astrophysical situations. Turbulent clumps evolve so 
that their column densities are equal to a critical value determined by the local FUV field, and 
typically have short photoevaporation timescales, ~ 10 4-5 years for a 1 M© clump in a typical 
star-forming region (r/ c o — 10, v — 10). Clumps with insufficient magnetic pressure support, 
and in strong FUV fields may be driven to collapse by the compressional effect of converging 
shock waves. We also estimate the rocket effect on photoevaporating clumps and find that 
it is significant only for the smallest clumps, with sizes much less than the extent of the 
PDR itself. Clumps that are confined by an interclump medium may either get completely 
photoevaporated, or may preserve a shielded core with a w r arm, dissociated, protective shell 
that absorbs the incident FUV flux. We compare our results with observations of some well- 
studied PDRs: the Orion Bar, M17SW, NGC 2023 and the Rosette Nebula. The data are 
consistent with both interpretations of clump origin, turbulence and pressure confinement, 
with a slight indication for favouring the turbulent model for clumps over pressure-confined 
clumps. 

Subject headings: ISMxlouds - ISM:dynamics - stars:early-type - stars: formation - 
ult,raviolet:stars 


1. Introduction 

Young massive OB stars significantly influence the structure, dynamics, chemistry and thermal 
balance of their associated molecular clouds through the impact of their ultraviolet photons. Their 
extreme ultraviolet photons (EUV, hu > 13.6 eV) ionize the gas immediately surrounding them, 
forming H II regions. Far ultraviolet photons (FUV, 6 eV < hu < 13.6 eV) dissociate molecular gas 
beyond the H II region, creating a photodissociation region or PDR. (Hollenbach & Tielens 1999). PDRs 
are ubiquitous: FUV photons dominate the gas heating in PDRs and affect the chemistry and physics 
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of gas over a large fraction of the volume and mass of molecular clouds. The study of interactions 
between FUV radiation and molecular cloud gas is therefore important in understanding molecular 
cloud evolution and the feedback between massive star formation and subsequent star formation in 
molecular clouds. 

Observations of PDRs compared with theoretical models probe the physical conditions in 
star-forming regions through dust and gas emission mainly in the infrared (IR) wavelength region of 
the spectrum. PDRs are the source of most of the IR emission in the Galaxy. Dust grains and large 
carbonaceous molecules such as polycyclic aromatic hydrocarbons (PAHs) absoib radiation from the 
stars, and re-radiate this energy flux in the infrared, producing a continuum with solid state and PAH 
spectral signatures. FUV photons incident on dust grains and PAHs also cause photoelectric ejection 
of electrons which then collisionally heat the gas, a process called the grain photoelectric heating 
mechanism (e.g., Watson 1972, Bakes & Tielens 1994, Weingartner & Draine 2001). The gas cools via 
emission in many IR lines. Emission from PDRs has been quite successfully modelled (see Hollenbach 
& Tielens 1999 and references therein). The models are usually parameterized by G 0 , the ratio of the 
FUV flux to the Habing (1968) FUV band flux of 1.6 x 10~ 3 erg cnrV 1 characteristic of the local 
interstellar radiation field, and n, the density of the gas. For advecting (non-stationarv) models, the 
flow velocity of the material through the PDR is also a parameter. The FUV fields in PDRs near OB 
star-forming regions such as the Orion Bar or Mlc SW are typically very high, (Go ~ 10 ), the PDR 

gas in these regions is inferred to be very dense with average densities ( n ) ~ 10 1 5 cm 3 (Tielens & 
Hollenbach 1985b). On larger sc ales and in more evolved regions such as the Rosette Nebula, the FUV 
fields and the average densities may be significantly lower by factors of 10 2 - 10 3 . 

Infrared, sub-millimeter and millimeter wavelength observations of PDRs indicate that gas in 
PDRs is inhomogeneous (van der Werf et al. 1993, Luhman et al. 1998), which is not surprising 
since molecular clouds themselves are observed to be highly inhomogeneous on a wide range of scales 
from tens of parsecs down to tenths of a parsec (e.g., Genzel 1991). There is both indirect and direct 
evidence for dumpiness in PDRs. Spatially extended fine structure line emission of neutral and 
singlv-ionized carbon near H II regions shows that the FUV penetrates deeper into the cloud than 
predicted by homogeneous models, suggesting that the gas is clumpy in nature (Meixner et al. 1992, 
Stutzki et al. 1988, Steiman-Cumeron et al. 1997). Further, strengths of high excitation lines of 
molecular species like CO J=14->13 are observed to be greater than that predicted to arise from gas 
with density inferred from other PDR species. It is thus inferred that in regions such as Ml 7 SW and 
the Orion Bar, the PDR gas consists of a “low” density (typically ~ 10 3 - 5 cm- 3 ) component responsible 
for the extended Cl and CII emission and high-density clumps (n ~ 10 6 cm ) which give rise to the 
high excitation CO lines from their surfaces or, more precisely, at a depth of A v ~1 from their surfaces 
(Burton et al. 1990, Meixner et al. 1992, Steiman-Cameron et al. 1997, Hogerheijde et al. 1995). 
Clumps are also observed directly by mapping PDRs at high spatial resolution in the 2 gm lines of H 2 
(e.g., Orion Bar, Luhman et al. 1998; Eagle Nebula, Allen et al. 1999). The small regions of enhanced 
emission are interpreted as indicating the presence of high density clumps (n ~ 10 6 7 cm ) with sizes 
< 0.02 pc, occupying a very small volume filling factor ~ 1 — 2%. How'ever, this interpretation is not 
very conclusive. Small differences in temperature within the observed regions can also explain the 
observed emission, as H 2 emission is highly sensitive to gas temperature (Marconi et al. 1998). 

Dense PDRs like the Orion Bar are spatially thin structures and if these PDRs are indeed clumpy, 
the clumps must have sizes smaller than the thickness of the w^arm PDR, ~ 0.03 pc for the Orion 
Bar. The presence of small clumps affects the infrared spectrum from PDRs by introducing a range 
of densities in the FUV-illuminated region and by introducing enhanced advection in clumps that are 



- 4 - 


photoevaporating. One of the goals of this paper is to investigate the conditions under which small, 
high-density clumps exist in PDRs. 

Small clumps might not exist in significant numbers if the intense FUV field photoevaporates them 
quickly. The incident FUV field heats up and pressurizes the surface layers of the clumps, causing 
the layers to expand. This mass loss may cause a complete photoevaporation of the dense clumps, 
and thus destroy them on short timescales. The mass loss timescales due to photoevaporation of an 
FUV-exposed clump can be estimated by assuming that the heated gas at the surface flows outwards 
at its sound speed, c PDH . The mass density at the base of this heated flow, p b , is lower than the mass 
density, p c , in the cold clump, as we will discuss below. (See Table 1 for a list of frequently used 3 
symbols in this paper.) The mass loss rate is given by dmjdt ~ 4?r r 2 c p b c PDR and since m c ~ p e r e , 
the photoevaporation timescale is t PE ~ rnj{dmjdt) ~ r c pJ{c PDR p b ). For a typical clump, and an 
FUV-heated surface at 1000K, this is approximately 3 x 10 3 (r c /0.01pc)(p c /p 6 ) years. Asymmetric 
mass loss from the surfaces of very small clumps can rocket them (e.g., Oort k Spitzer 1955) out of 
the PDR back into the cloud. Clumps that survive the rocket effect and remain in the PDR are likely 
to lose significant fractions of their mass on short timescales of order J 10 \ ears, depending on 

the ratio pj Pb- The clump evaporation timescale can be larger or shorter than its residence time in 
the FUV-illuminated surface region of a CMC. From a frame of reference moving with the advancing 
ionization front, interclump material and the clumps in the GMC are advected into the PDR and 
ultimately flow into the H II region. Large clumps survive and they can penetrate the H II region and 
affect its evolution and structure, becoming Evaporating Gaseous Globules or “EGGs (Bertoldi k 
McKee 1990). 

In the present paper, we focus on the heating of an individual clump in a PDR by FUV photons 
from a nearby O or B star. We study a range of FUV fluxes incident on a range of clump sizes and 
determine the photoevaporative lifetimes of the clumps and the evolution of their structure. Our 
investigation is analogous to the study by Bertoldi (1989) and Bertoldi and McKee (1990) on the 
effects of EUV heating of clumps in H II regions, and a generalization of the studies by Johnstone et al. 
(1998) and Storzer k Hollenbach (1999) of the evaporation of small clumps and protoplanetary disks 
( “proplyds” , O’Dell et al. 1993} by EUV and FUV photons. 

In a subsequent paper (Gorti k Hollenbach, in preparation) we will discuss the cumulative 
effect of O and B stars on Giant Molecular Clouds (GMCs), and the relative importance of EUV 
and FUV photons in dissociating and destroying GMCs. Ultimately, we are interested in the role 
of FUV radiation in regulating star-formation in GMCs. Destruction through photodissociation and 
photoevaporation of potentially star-forming clumps/cores in GMCs reduces the star-forming efficiency 
of GMCs, and the FUV field may thereby play an important role in regulating or limiting star 
formation. Photoevaporation of the outer envelope of a collapsing isothermal sphere by FUV radiation 
could rapidly deplete the material in the envelope, and may limit the final mass of the star formed. On 
the other hand, compression of non-collapsing clumps by shock waves driven by the warm suiface gas 
could possibly drive the inner cores to instability and gravitational collapse, triggering star formation. 

In §2 of this paper, we discuss the present understanding of GMC structure and existing models. 
In §3, we describe the evolution of an H II region and its surrounding PDR as they propagate into the 
GMC. Timescales relevant for the FUV heating of a clump in various contexts are briefly discussed 
in §4. In §5 w r e present a simple analytical model which provides an understanding of the underlying 
physics of FUV-induced clump evolution. We consider various possible evolutionary scenarios and 
their dependences on clump parameters and the strength of the FUV field. For the more general case, 
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a 1-D numerical hydrodynamics code has been developed and these results are described (§6). In §7, 
we estimate the significance of the rocket effect due to asymmetric mass loss and its effect on clump 
lifetimes in the PDR, and in §8 we discuss the implications of our results for observations of PDRs m 
star-forming regions. Application to observed PDRs is made in §9. We conclude with a summary of 
the present investigation (§10) 


2. Clumps in Giant Molecular Clouds 

Giant Molecular Clouds are observed to have a hierarchical substructure that is highly 
inhomogeneous, with density enhancements such as cores, clumps and filaments on all observable 
scales from ~ 100 pc down to ~ 0.1 pc. We use the term “ massive cores” to denote large scale density 
enhancements (> 1 pc in size) and “clumps” to denote smaller scale structures (~ few tenths of pc 
or less in size^GMCs and their more massive constituent cores are strongly self-gravitating, and 
are supported against collapse by turbulence and magnetic fields. The internal velocity dispersions of 
the substructures in molecular clouds are found to scale with their sizes, with decreasing non-thermal 
support on the smallest scales (Goodman et al. 1998). The column density through a GMC or a 
typical GMC substructure is found to be approximately a constant N ~ 10 22 cm , although significant, 
variations (N ~ 10 21 " 23 cm- 2 ) about this typical column may occur. (For a recent review on GMCs, 

see McKee 1999). 

The smaller substructures or clumps are not all gravitationally bound and the exact physical 
nature of these clumps remains uncertain (Williams et al. 2000). They have been interpreted as 
temporary fluctuations in density caused by supersonic turbulence within the cloud (e.g., Falgarone 
& Phillips 1990, Scalo 1990). This view is supported by their apparent fractal nature, which seems 
to show the same structure on smaller scales down to a few tenths of a parsec (Elmegreen 1999), and 
their observed supersonic linewidths (Blitz & Stark 1986). In this picture, the gravitationally unbound 
clumps are transient objects formed by converging supersonic turbulent flows. The non-gravitating 
clumps are therefore being constantly formed and destroyed, and have short lifetimes of the order 
of their sound crossing time, t c « r c0 /c c , where c c is the clump sound speed, and r c0 is the initial 
clump radius. The effect of the sudden turn-on of an FUV field on such clumps is to drive a shock 
into the clump, producing a smaller, denser clump with a warm evaporative outflow. Surprisingly, the 
evaporative timescale is now somewhat longer than t c , as we will show. Assuming that the FUV flux 
does not strongly affect the turbulence so that the formation timescale remains the same, the presence 
of an FUV field will not greatly change the non-gravitationally bound clump abundance, but it, will 
compress them and alter their thermal structure and dynamical evolution. 

Another interpretation is that small bound clumps as well as small unbound dumps represent 
stable physical entities confined by interclump pressure. Observations through molecular line spatial 
maps and velocity channel maps or position-velocity diagrams indicate the presence of dense coherent 


iNote that our terminology of cloud substructure is the reverse of what is sometimes used in ^ the 
literature. We call the large scale structures “ massive cores” and the small scale structures clumps to 
be consistent with observational work on dense regions in PDRs, wherein the density enhancements are 
usually called “clumps” . Our notation is also consistent with the notion of “hot cores” , where massive 

stars form. 
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structures. Indirect evidence suggests the presence of a low density (~ 1 - 2 orders of magnitude 
less than the clumps) “interclump medium” (ICM) (see Williams et al. 1995 and references therein). 
The ICM exhibits high velocity motions (Au « 10 km s _1 ) and therefore may have significant 
pressure, which could be of thermal, magnetic or turbulent origin. Clumps in this picture are confined 
by the ICM, and some of them may be self-gravitating (Maloney 1988, Bertoldi k McKee 1992). 
Pressure-confined clumps live much longer and presumably form more slowly than the gravitationally 
unbound clumps of the turbulent model. Therefore, FUV heating of these clumps and their destruction 
through photoevaporation will have a much more pronounced effect on lowering the clump abundance 
in this interpretation. 


3. The Overall Evolution Of A Massive Star Forming Region 

Massive stars form in the central regions of clumpy, massive, cores (M > 1OOOM 0 ), of radius ~ 1 
pc inside GMCs (see review by Evans 1999). These massive cores contain clumps whose densities 
may reach n c ~ 10 5 ~ 7 cm -3 , or even higher (Plume et al 1997), and the clumps are surrounded by 
interclump gas with lower, but uncertain densities, possibly n IC M ~ 10 3-5 cm -3 . W hen a massive star 
forms, the EUV photons from the star immediately ionize the nearby interclump gas, and form an H II 
region. The clumps in the H II region are also exposed to the ionizing radiation and evaporate (Bertoldi 
1989, Bertoldi and McKee 1990). The FUV photons penetrate further through the interclump gas 
beyond the interclump ionization front (IF) and dissociate the molecular gas, forming a dissociation 
front (DF) and a PDR. 

In the earliest stages of evolution when EUV photons emerge from the forming O/B star, the 
IF moves as an R-type front. However, the IF quickly stalls and becomes D-type, at which point 
the high pressure H II region drives a ~ 10 kms -1 shock into the neutral interclump gas beyond 
the IF. It is instructive to obtain a sense of sizescale at the R/D type transition. Assuming the 
massive star emits fa Lyman continuum photons per second, the Stromgren radius at this juncture is 
r s « 0.1<^V 4 2/3 pc, where fa 9 = fa/10 49 ^ 1 and the electron number density n e4 = n e J 10 4 cm~ 3 . If 
the central star has an FUV luminosity of L FUV = 1O 5 L 0 (typical of stars with fa 9 « 1), the FUV 
flux incident on the surrounding PDR is characterized by G 0 * 10 5 . The thickness of the PDR is 
approximately 10 22 cm- 2 /n /C M or 0.3 nJ^ M4 pc. The DF initially advances ahead of the shock, but 
as photodissociation approaches equilibrium with H 2 formation, the DF slows down. Eventually, the 
shock moves ahead of the DF, and the entire PDR is contained in the post-shock shell. The shock 
now impacts cold molecular interclump gas. Hill k Hollenbach (1978) showed that at this stage the 
shock velocity has slowed to < 3km s" 1 , and that, if n ICM > 10 4 cm" 3 , the shock completely stalls in 
t < 10 5 years. Therefore, there is an initial stage {t < 10 5 years, for dense ICM) in the evolution of 
the clumpy PDRs surrounding H II regions where the clumps are both shocked and exposed to the 
photoevaporating effects of the FUV radiation. This can be followed by a stage where the shock has 
stalled and unshocked clumps move into the FUV-illuminated zone. 

Two configurations are possible after the shock stalls. If the H II region is still completely 
embedded, it becomes stationary with its pressure matched by the molecular cloud pressure and the 
ionizing photons completely absorbed by recombining electron/proton pairs in the H II region. Much 
more likely is the situation where the H II region breaks through the surface of the cloud and becomes 
a “blister” H II region. In this case, the ionized gas can expand away from the cloud, and into the 
interstellar medium. The IF then slowly eats its way into the PDR and the DF similarly advances into 
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the molecular cloud. The velocity of the IF advancing into the PDR is of order v f ~ 1 kms , just 
enough so that the flux of particles through the IF can balance the heavily attenuated flux of EU v 
photons reaching the IF (see Whitworth 1979, Bertoldi k Draine 1996, Storzer k Hollenbach 1998). 
In such a case, the clumps entering the PDR may not have been previously shocked and compressed. 
Further, in the turbulent GMC model, clumps form and dissipate continually and those formed in the 
PDR after the passage of the shock are not affected by it. In both the PDR around the Trapezium 
stars in Orion and in the M17 SW PDR (Meixner et al. 1992), there is no evidence for velocity shifts 
between PDR gas and the ambient cloud gas, or, in other words, no evidence for a shock at the present 
time. In this paper, we focus only on the photoevaporation of clumps, whether or not they have been 
shocked, and do not follow the potential shredding or flattening of a clump by the passage of a shock 
wave (see, e.g. Klein et al. 1994). 

The thickness or column depth of a PDR has been defined by Tielens k Hollenbach (1985) to 
include all the gas and dust where FUV photons play a significant role in the gas heating or chemistry. 
By this definition, in regions like Orion or M17, the PDR extends to columns of - 10 cm and 
includes gas which is cool and molecular (H 2 , CO), but where FUV still photodissociates 0 2 and H 2 0. 
W 7 e are interested in the surface column N 0 of the PDR which is heated to the highest temperatures 
(T maI ~ 100 - 3000K, depending on G 0 and n). For G Q /n > 4 x 10~ 2 cm 3 , N o «l-3xl0 cm 
(equivalent to a dust FUV optical depth of order unity) and is nearly equal to the column where 
atomic hydrogen converts to H 2 in stationary PDRs (Tielens k Hollenbach 1985). We shall call this 
region the “warm PDR” hereafter to distinguish it from the entire PDR. The thickness of the warm 
PDR, No/n/cAr, is of order 0.03 - 0.1nj£ M4 pc. The crossing timescale for the IF through the warm 
PDR is t PDR « 0.3 - 1 x 10^ n}c M4 vJs years, where v fb = v f /(l km s *). Gas and dust (and clumps) 
from the opaque molecular cloud interior are therefore advected on these timescales from the shielded 
cloud, through the warm PDR. and across the IF into the H II region. Figure 1 shows a schematic 
diagram of FUV radiation from nearby O/B stars impinging on massive cores (large sizescales, e.g., 
Rosette) and clumps (small scales, e.g., Orion). 


4. Impulsive heating versus slow heating 

FUV-illuminated clumps are heated at their surfaces on timescales t FU v that could be either fast, 
or slow compared to the internal sound crossing timescale w'ithin the clump, t r = r c o jc c . Impulsive 
heating implies that the FUV heating of the clump surface layer, characterized by column N 0 , occurs 
with t FU v < t c . In this case, the outer region is heated more quickly than the entire clump can 
respond, and attains a much higher pressure than the cold central region. A shock propagates inward 
and a heated photoevaporative flow expands off the surface. 

Several situations in clumpy molecular clouds may lead to t F yy t c . For a clump with r c o ~ 10 
cm and c c ~ 0.3 km s" 1 (T ~ 10K), t c ~ 10 4 5 years. The OB star may “turn-on” on a timescale <10 
years. Even after this initial flash of FUV flux, the clumps may shadow each other so that t F uv is given 
by the time for a clump to move out of shadow, which is of order r c0 /v c , where v c is the clump velocity. 
In a turbulent, supersonic medium v c may be larger than c c so that t FU v < t c . Similarly, if clumps 
are formed in turbulent gas by supersonic converging flows, the FUV turn-on timescale is roughly the 
time for the clump to form, which again is t F yy ~ 'fco/ v ti where v t is the supersonic turbulent speed. 
Therefore, impulsive heating is likely the most appropriate approximation for the turbulent model of 
clumps in molecular gas. 
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Let us consider the other situation of long-lived clumps confined by the ICM. Here, slow heating 
is generally the best approximation, but in certain conditions, impulsive heating may be more 
appropriate. A common situation after the turn-on of the OB star may be the advancement of the 
ionization front and PDR front into the opaque molecular cloud, which occurs at speeds v f ~ 1 km 
s _1 . In this situation t F pv is the time for a PDR clump to move one FUV optical depth closer to the 
ionization front, or tpuv ~ tpDH = Xpdr! v Ji where X .pdr — Xq/picm is the thickness of the heated 
PDR zone in the ICM of density ti, C m- The criterion for impulsive heating t FUV < t c , is more easily 
met by larger clumps that have longer sound crossing timescales, as t FUV is the same for all clumps. 
This gives us a limiting size for clumps likely to be impulsively heated, 

r c o > Xpor — ~ (0.3 — 1)Xpdr, (1) 

v f 

which is of the order of the thickness of the interclump PDR region itself. Clumps with radii smaller 
than this are slowly heated as they emerge from the opaque molecular cloud into the FUV-irradiated 
regions. Since we only consider clumps within PDRs that have r c o < X 'pdr, pressure-confined clumps 
are generally slowly heated in our analysis. 

Clumps that undergo slow heating, t F uv adjust quasi-statically to the changing FU\ 

flux. Long-lived, small, pressure-confined clumps advecting into FUV-irradiated regions provide a 
prime example of this slow evolution. Such clumps develop a. warm PDR surface which slowly gets 
warmer, thicker and less dense. The clump and the heated PDR surface evolve in near pressure 
equilibrium, with the density of the surface decreasing with the rise in its temperature, such that 
ppdrTpdr = n c T c = pressure of the ICM (assuming thermal pressure dominates on these small scales). 
At any given time, the evolutionary state of a slowly heated clump is the same as the final state of an 
impulsively heated clump exposed to the same local FUV flux. However, a gradually heated clump 
constantly adjusts its PDR surface to maintain pressure equilibrium as the PDR slowly heats. In 
contrast, an impulsively heated clump undergoes shocks and supersonic expanding flows in its attempt 
to attain equilibrium. (We note that there is no final equilibrium state if the interclump pressure 
is zero). Rapid photoevaporative mass loss can cause an acceleration due to a rocket effect in very 
small clumps and move them back into the shielded cloud. Pressure-confined small clumps, on the 
other hand, are not rocketed as they evolve quasi-staticallv. In a given FUV field and at a given ICM 
pressure, the final equilibrium state of a small clump is either a completely photodissociated, warm 
(T PD r), and expanded (but with P = P, C m ) region; or it may consist of a shielded cold molecular 
core with the same temperature, density and pressure as the original clump, but with a heated and 
more diffuse protective surface layer of PDR material. In §5.3 we provide analytic solutions for these 
equilibrium structures. However, in §5 and §7 we primarily focus on the interesting physical processes 
which occur for clumps that are impulsively heated. 


5. Analytical models for FUV photoevaporating clumps 

We first provide approximate analytic solutions for the time evolution of a clump suddenly exposed 
to FUV radiation, using a few simplifying assumptions. As described in §2, the nature of the clumps 
themselves is not well understood, and we consider two simple analytic models that qualitatively 
correspond to the two main interpretations. Turbulent unbound clumps are modelled as constant 
density structures in vacuum. A clump in vacuum disperses on a sound crossing timescale, in the 
absence of FUV illumination. Pressure-confined clumps are also treated as being of constant density, 



but have a surrounding ICM with equal pressure. In our models, the ICM only serves to confine the 
clumps and does not itself evolve or otherwise affect the clump evolution. We seek solutions for the 
time evolution of the mass, size and mass loss rate of a photoevaporating clump. 


5.1. Initial conditions and assumptions 

Clumps are assumed to be dense, small spheres of gas, supported by thermal, turbulent and 
magnetic pressures. We assume that the magnetic field B scales with a constant power of the density, 
so that the magnetic pressure is 

P B = B 2 /8nxn\ (2) 

The initial ratio of turbulent and magnetic pressures, P NT and P B , are set by two dimensionless 
parameters, 

n = P nt / Pt\ ft = Pb/ Pt (3) 

where Pp is the initial thermal gas pressure in the clump. A clump may have significant turbulent 
support and observations of clumps in star-forming regions (Jijina et al. 1999) indicate values of 
q; ranging from ~ 0 in cold, dark clouds to ^ 2 in regions of massive star formation. There is 
observational evidence to suggest that turbulent support in clumps decreases on the smaller scales of 
star-forming clumps (r c0 < 0.1 pc), where a < 1 (Goodman et al. 1998). Measurements of magnetic 
fields and hence ft and 7 are difficult, but present observational data suggest a wide range of values for 
from ~ 0 to a few (Crutcher 1999). The value of 7 (~ 1 - 2) describes how the B field responds to 
a relatively sudden change in gas density. For a frozen, uniform magnetic field threading a spherical 
clump, B oc n 2 P or 7 = 4/3. For our standard case we adopt values of a — 1 ,8= 1, and 7 = 4/3. 
We later discuss the sensitivity of our results to these choices for a, ft and 7. As discussed in §4, for 
clumps confined by an ICM and adverting into FUV-illuminated interclump zones, the assumption of 
tFUv/tc < 1 is equivalent to r c 0 > (0.3 — l)X pdr- Therefore, in this case we treat only clumps above a 
certain minimum size, which is of order the size of the thickness of the heated PDR of the ICM. We do 
not explicitly include self- gravity in our analytical model, and only qualitatively estimate the effects of 
gravity on clump evolution. We treat gravity more quantitatively in the numerical models. 

The FUV radiation field as seen by the clump is spherically asymmetric, and stronger at the 
surface facing the source. For simplicity, however, we assume a spherically symmetric FUV field when 
calculating the internal dynamical evolution and the thermal structure of the evaporating clump. 
Typical albedos of interstellar dust are ~ 0.5 and some of this scattered radiation is directed backwards 
(dependent on the phase factor c/, see Henvey & Greenstein 1941), resulting in a backscattered FUV 
flux on the shielded clump surface about 0.1 - 0.2 times that on the directly illuminated surface (e.g., 
Desert et al. 1990, Hurwitz et al. 1991). Although the posterior intensity is diminished compared to ^ 
that on the front surface, the PDR surface sound speed is not a sensitive function of the incident FUV 
flux G 0 (Kaufman et al. 1999). For field strengths G 0 - 10 3 ” 5 near typical H II regions and densities 
n > 10 3-6 cm” 3 characteristic of clumps there, even a back flux 10% of that on the clump surface 
facing the source would heat the gas to nearly the same sound speed (cp BB — 3 km s 1 on the source 
side, cpoR — 2 km s” 1 on the shadowed side). Therefore, the assumption of spherical symmetry is as 
good as ignoring the flux of FUV photons on the backside when calculating the thermal structure and 
internal dynamics. However, the asymmetry of the flows from the front and back surface will cause a 
rocket effect on the clumps and induce motion away from the source, as first recognized by Oort & 
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Spitzer (1955). We discuss momentum transfer due to mass loss of the clumps in §7, and estimate sizes 
of clumps that gain velocities high enough to move out of the PDR and back into the molecular cloud. 

In our analytical models (but not our numerical models) we assume that the transition from the 
warm outer PDR layer to the cold molecular core of a clump is very thin, so that the clump can be 
modelled as an isothermal core with sound speed c c , surrounded by a PDR envelope of sound speed 
c PD R. The gas in the core always remains isothermal in our analysis, even when shock-compressed to 
high densities. For our pressure- confined clumps, it is assumed that the ICM is of low density and does 
not absorb any FUV photons. In order to obtain analytic solutions, we do not fully treat transient 
phenomena and emphasise quasi-steady state aspects of the flow. 


5.2. Analytic model for unconfined clumps 

We begin by considering the evolution of a simple configuration: a cloud consisting of clumps that 
are spheres of radius r c0 and constant initial density n c0 , immersed in a vacuum. We refer to such 
clumps hereafter as “turbulent clumps” . These are dynamic structures that in the absence of FUV 
heating expand at roughly their sound speed, c c . When such a clump is exposed to an FUV field, the 
photon flux is attenuated by dust and instantaneously heats an outer column of N 0 ~ 2 x 10 21 cm 2 
(thickness ($o = -Vo/^co) to relatively high temperatures and sound speeds (for example, near H II 
regions with Go ~ 10 4- ® and ti c q ~ 10 5 cm 3 , T ~ 1000K and cpop ~ 3 km s ). 

The evolution of the clump is determined by its initial radial column density from the centre 
outwards, fico^co-i and the ratio u of sound speeds in the FUV -heated region and the cold clump 
material, 

v = Cpdr/c-c ( 4 ) 

The parameter v can be thought of as a measure of the strength of the FUV field incident upon the 
cold clump. We introduce the dimensionless parameter r] c o, 

VcO = n c or c o/N 0 - r c0 /6 0 , (5) 

for the initial column density to the centre of the clump. Note that t) c q is a measure of the mean 
column density through a clump. 

For turbulent clumps with 7 ] c q < 1, the initial clump column density is less than iVo, and the entire 
clump is immediately heated and photodissociated by the FUV flux. Effectively, therefore, the FUV 
flux accelerates the expansion of the cold clump by heating it throughout, decreasing the expansion 
timescale by a factor given by the ratio of sound speeds at the two temperatuies, c c j cpdr, or 1//a 
For clumps with t] c q > 1, the FUV flux heats an outer shell of gas to a higher temperature, thereby 
increasing its pressure. If the pressure of the outer PDR shell is sufficiently high, it drives a shock to 
the centre of the clump rapidly compressing it, and the compressed clump proceeds to evaporate on 
a somewhat longer timescale than the expansion timescale in the absence of an FUV field. On the 
other hand, if the pressure in the heated outer PDR shell is low, the shock does not make it to the 
centre of the clump. The shock dissipates, followed by an expansion of the clump, until finally the 
photoevaporative flow halts the expansion and proceeds to shrink and compress the gas. There are 
thus two distinct evolutionary scenarios for photoevaporating turbulent clumps with t] c q > 1. 

For brevity in the main text, the details of our analysis and the derivations of the results are 
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presented in the appendices to the paper. We now discuss our solutions and describe the physics of 
FUV-heated clump evolution and refer the interested reader to the appendices for the full analysis. 

Clumps exposed to a given FUV field, measured by the parameter is, are either compressed b^ 
shocks or expand to adjust their column density to a critical value, Vcrit- This critical column density, 
7 j crit can be derived from conservation of mass and from the condition of pressure equilibrium at the 
surface of the cold clump gas (see Appendix A for details), and is defined by the relation 

2 ^ +a l + 0(—L—Y = i+ 0 + / ? (6) 

Went “1 V Wcrit 1 / 

Recall that a and are the initial ratios of turbulent to thermal pressure and magnetic to thermal 
pressure, respectively, and that the magnetic pressure scales as n? . For typical values of these 
parameters, the second term on the LHS is negligible and equation (6) can be used to define a 
photoevaporation parameter, A. as 

. _ (1 + ft 4- P)(jlco ~ i) (7) 

2(2is 2 + ft) 

For clumps with initial column densities r ) c 0 = rj cril , the photoevaporation parameter A = 1. Figure 2 
shows an t ] — is parameter plot with equation (6) depicted for the standard case with ft F/5 U 
and 7 = 4/3, and for two possible extremes of these parameters. The result is not very sensitive to 
the choice of the parameter 7, and Vcrit only depends on the sum of ft and /3. for the standard case 
and with t] c0 » 1 and is » 1 which is typical of many PDRs, equation (6) can be approximated 
as Vent ~ 4n 2 /3. The evolution of FUV-heated clumps is determined by where they initially lie 
in the 77 — is parameter space. Clumps with an initial column density tj c q < Vcrit evolve through a 
shock-compression phase, until their dimensionless column density increases to Vcrit (see Appendix B). 
Clumps with higher initial column densities expand until their column density decreases to rj crit (see 
Appendix C). The evolution of clumps in two representative cases is qualitatively shown in Figure 2 
through their trajectories in 7/ -- is space. Once the clump column density attains its critical value T] cr n, 
further evolution is similar in both cases, and the clump column density at any later time remains 
constant at this critical value. We note that this would imply that observed clumps in a given PDR 
should typically have the same column density, the critical value r) crit given by the ambient FUV field. 
This is valid for clumps photoevaporating into a vacuum and only violated in the transient early stages 
which are typically very short, of order t. « r c0 /c PDR . From equation (6), clumps at 10K heated on 
their surfaces to 100 K (or is ~ 3 ) should have columns VcritNo ~ 3 x 10 22 cm 2 , and for T ~ 1000/v 
or is ~ 10, Tjcrit No ~ 3 x 10 23 cm -2 . As photoevaporation proceeds, and clumps lose mass, they shrink 
in size and their average densily increases, so that the column through the clump remains the same, 
n c (t)r c (t) = r) crit N 0 . 


Shock compression In this case, the initial column density is larger than No, but smaller than 
Vent No- The FUV can only penetrate through an outer shell of thickness S 0 = No/n c0 < r c0 . This shell 
is heated to a higher temperature and the corresponding increase in pressure causes the outer shell to 
expand at cpp R radially outwards at the edge and inwards into the clump at the inner boundary. Figure 
3 shows a schematic diagram of the evolution of a shock-compressed clump. As the shell expands, the 
column density N 0 through it is maintained by further penetration of FUV photons, and heating of 
more cold clump material. The pressure at the base of the PDR shell is initially higher than the cold 
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clump gas. This pressure difference drives a. shock wave into the clump, which rapidly propagates to 
the centre Behind the shock the cold clump gas is compressed to a pressure approximately equal to 
the pressure at the base of the PDR shell. The shock travels inward at roughly c PDR and the entire 
clump gets compressed in a time t s ~ r c (0)/c . PDR , where r c (0) = r c0 - S 0 is the initial radius of t e 
cold clump when the FUV flux is turned on (Fig. 3b). The shock reaches the centre and compresses 
the clump to a radius r c {t s ) = r.,. The shock-compressed gas forms a high-density core, such that the 
pressures in the core and at the base of the PDR flow are approximately equal. Further mass oss from 
the clump is from the surface of this core. The time evolution of the mass and radius of the clump tor 
t > t s are given by the following equations for our standard case, where a — 1, /? — 1, and 7 ~ 4/ 

For the general solution, refer to Appendix B. 


r c {t > t s ) 


r$ 

r c{0) 


5/4 


10 ^ 9/4 ( 1 

3(l7cO - l) 2 V*C 


{rjco ~ 

W ) ) 


4/5 


r c (0) 


( 8 ) 


drric _ f mco\ f ^c(0\ (9) 

dt ~ \ t c ) \ r c o ) rjco 

m c (l>t,) = m c ( 0)(^y) ( 10) 

Here q = (A/3) 1 / 3 , r s « qr c { 0)(1 - 3(1 + q)/{VcO ~ 1)) 4/9 and t > t s = r c (0)/c PDR . The PDR/core 
interface propagates inward as the warm gas evaporates and, finally, the entire core is transformed 
to heated, photodissociated, expanding PDR material. Using equations (8) and (10), Figure 4 shows 
the clump radius and mass as a function of time for a case where the initial column through the 
clump is given by Vco = 10, and v = 5 represents the incident FUV flux. This is a typical column 
density for clumps in molecular clouds, (see §2) and for a clump of mass m c , corresponds to a density 
n = 3 x 10 6 (M 0 /m c ) 1 / 2 cm -3 . In this example of v = 5 , the FUV field might heat, cold clump gas 
with initial temperature of ~30 K to about ~750 K, or ~ 10K gas to about 2o0 K. Also shown m 
Figure 4 are the results of a more detailed hydrodynamics code simulation (described in the next 
section) for the same parameters. There is very good agreement between the two results, in spite of the 
various simplifying assumptions made in arriving at the analytical solutions. The numerical solution 
for r c (t) oscillates with time due to the dynamics of the shock compression which we have neglected 
in the above analysis. Our analytic solution assumes a quasi-steadv state, whereas the numerical code 
shows the effect, of overshooting the steady state solution and rebounding. 


The clump evolution goes through two distinct phases, which is apparent from the plot of clump 
radius with time (Fig.[4a]). The dump initially shrinks rapidly with time, as the shock compresses 
the clump in a time t s . The average velocity v b with which the radius decreases is 0.65 c PDR for the 
chosen parameters. After the entire clump is compressed by the shock, it decreases in size more slowly, 
now due entirely to the mass loss from its surface. As discussed earlier, the column density in the 
compressed clump remains constant, and the clump gets denser as its size decreases. In Figure 4b, 
the mass of the clump is also seen to decrease with time as the entire clump gets photoevaporated in 
about a sound crossing time; again, the analytical results closely follow the result of the numerical 
hydrodynamical code. The cold clump mass at t = 0 is less than the initial clump mass because it 
excludes the mass contained in the outer column N 0 which is instantaneously heated by the FUV flux 
to the PDR temperature. 

Using the results of the above analvsis, we can obtain simple estimates of photoevaporation 
timescales for a clump exposed to incident FUV radiation. We define the photoevaporation timescale, 



- 13 - 


tpE, as the time for the radius of the clump to shrink to zero. The time can be easily determined 
from setting the LHS of equation (8) to zero, 

tpE ~ 0.5^ 3 i/ -1/3 t c (11) 


or 


tpE ~ 10 


A( 


n, o 

1 10 5 cm -3 


2/3 


0 


O.Olpc 


5,/3 / 0.3kms 1 \ ^ ( 3kms 1 \ 


) e 


CPDR 


years 


(12) 


Larger clumps are longer-lived compared to smaller clumps in a given PDR environment. 
Photoevaporation timescales of turbulent clumps in PDRs in typical star-forming regions are thus of 
the order of a few clump sound crossing times (see §7). Paradoxically, the sudden turn-on of FUV 
radiation on a clump in a vacuum can increase its lifetime, even though the FUV heats the surface 
and increases the flow speed from the surface! For turbulent clumps with lifetimes ~ t c in the absence 
of an FUV field, exposure to the FUV field results in a “pressure-confined” period where the clump is 
compressed due to pressure of the heated surface PDR layer. This compression reduces the area of the 
photoevaporating surface and can extend the lifetimes of these clumps to a few t c . 


Collapse of clumps driven by shock compression The FUV driven shock wave that compresses 
a clump to very high densities, may render it gravitationally unstable to collapse. Thus star formation 
may be triggered in a previously stable clump. A simple estimate of parameters leading to clump 
collapse can be made without explicitly including self-gravity in the equations. We use our solutions 
for the radius of the shock-compressed core (Appendix B) and compare this with its Jeans length, 
rj. If r s < rj, the clump is stable, otherwise it undergoes gravitational collapse. It should be noted 
here that a clump supported by magnetic pressure with 7 > 4/3 will always be stable to collapse 
regardless of the external pressure, as the magnetic pressure increases faster during compression than 
the gravitational energy density. Our discussion here is thus applicable for clumps with insignificant 
magnetic fields, or with a magnetic equation of state where 7 < 4/3. We now solve for the collapse 
criterion as a function of initial clump parameters and the FUV field strength as measured by u , for an 
illustrative case where a = l,/3 = 1, and 7=1. The radius of the shock-compressed core is given by 


r s ~ r c (0) 


/ 3fa„- 1 ) \ 1/2 

U^ + uy 


The Jeans length of the compressed core is given by 


rj 
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For collapse r s > rj which yields 


VcO > 1 + T 


/ 7rn c0 c 2 \ 

\Gm„N*) 


2/3 

(r/ 2 + I)- 1 / 3 


(13) 


(14) 


(15) 


The factor n c0 c 2 in equation (15) is the initial thermal pressure in the clump, and for clumps in PDRs 
is typically of the order of 10 6-7 cm -3 K. For n c oC 2 = 10 6 cm -3 K, equation (15) becomes 

49 


VcO > 1 T 


(l/ 2 + l) 1 / 3 


(16) 
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Figure 5 shows the collapse criterion(Eq.[16]) on the r)-v parameter plot. Clumps with initial column 
densities greater than ~ 40iV o and with initial thermal pressures of 10 cm ‘ K cannot support 
themselves against gravity and are unstable to collapse even in the absence of an FUV field_ is 
region of instability is to the right of the vertical dotted line in the plot. Clumps to the left of this me 
are initially stable; but sufficiently strong FUV fields may drive shocks that compress them and trigger 
collapse. The FUV fields near OB stars may thus trigger star formation in previously stable clumps, 
and increase the star formation rate in clouds, in cases where 7 < 4/3. 


Clumps with an initial expansion phase Clumps with initial column densities r] c0 > T] CTl t a& 
given by equation (6) develop very thin (So « r c0 ) PDR shells when they are first exposed to FUV 
radiation. This very thin shell is initially at a high pressure and drives a shock into the clump on the 
inside edge of the shell while the outside edge of the shell expands at roughly c PDR just as is the case 
for shock-compressed clumps. However, the shock stalls before it reaches the clump centre as the PDR 
pressure rapidly declines due to the shell expansion. The PDR pressure and clump pressure become 
equal before the shell thickness is comparable to the clump radius. We summarize below the detailed 
evolution of initially expanding clumps; details appear in Appendix C. 

Because of the high PDR sound speed and relatively small thickness, the PDR pressure initially 
drops faster than the pressure in the cold clump. The pressure in the cold clump eventually becomes 
higher than that in the surrounding expanding PDR shell, and the cold clump gas expands at its sound 
speed c c (as it would in the absence of FUV radiation since the interclump medium is assumed here 
to be a vacuum). This expansion takes place until a time t e at which the density (and pressure) in 
the clump gas drops to that of t.he outer PDR layer. For t > t e , FUV photons begin penetrating into 
the formerly shielded cold clump gas. The warm PDR gas now confines the clump gas, and further 
expansion of the cold clump is halted. The cold clump proceeds to lose mass gradually and shrinks due 
to photoevaporation, with an evolution similar to the final solution for shock-compressed clumps. The 
radius and mass of the clump av times t > t e and for o; = 1, 8 < 1, 7 < 4 / 3 are S Iven b >' 


r e (t > t e ) = r c (t e ) 
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where r c (t e ) = r c (0) + c c t e , 
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(18) 


and 


, ( 2^ 2 + 1 ) 

rn c (t > t e ) = rn c o 

VcO 



(19) 


Figure 6 shows the change in radius and mass of a clump with initial parameters r) c0 100, v 3. 
The clump expands out initially and loses a small fraction of its mass during this phase. After 
its density has dropped, so that there is pressure equilibrium between the cold clump gas and the 
heated layer immediately surrounding it, expansion is halted. As photoevaporation causes mass loss 
off its surface, the clump shrinks in size, and is completely heated and photodissociated in about 
4 crossing timescales. The numerical results are also overlaid on the plot, and there is reasonable 
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agreement. The photoevaporation timescale for initially expanding clumps can be defined analogously 
to shock-compressed clumps as '.see Appendix C, Eq.(C7)) 
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Therefore in this regime of low u, the photoevaporation timescale t PP is of order several sound crossing 
timescales, and is proportional to u. 


5.3. Analytical model for clumps pressure-confined by an interclump medium 

For a cloud model in which clumps are structures in pressure equilibrium with an ICM, the 
evolution of FUV-heated clumps is somewhat different. The heated shell or PDR can now no longer 
expand indefinitely, but is eventually confined by the pressure of the ICM. The evolution of the clump 
is again determined by ry c0 and v. Clumps with an initial column density r? c0 < N 0 , are completely 
heated and photodissociated instantly. Such clumps expand until the PDR gas reaches the interclump 
pressure. Clumps with higher initial column densities again are only heated on their surfaces to a 
column 7V 0 , and their further evolution depends on the turn-on time of the FUV field t FUV , relative to 

fe- 
lt the clumps are heated impulsively, t FUV <C t e , pressure-confined clumps initially evolve similarly 
to the turbulent clumps discussed earlier. Clumps with r/ c0 < T] crit are thus shock-compressed, and mass 
flows are set up which photoevaporate the clump. However, the outflow runs up against the interclump 
pressure and eventually reaches pressure equilibrium with its surroundings. During this time, the 
clump may either be completely heated through and get transformed into a sphere of warm PDR 
gas, or may only be partially photoevaporated. The final configuration for partially photoevaporated 
clumps consists” of a small remnant cold clump, and an extended warm PDR envelope surrounding 
it and protecting it from the FUV flux. The clump, PDR envelope and interclump medium are in 
pressure equilibrium. Clumps with r] c0 > i] cr it are not much affected by the FUV field. There is no 
expansion of the clump due to the presence of the confining ICM, and the initially heated thin PDR 
shell expands to form an expanded, but still thin, protective PDR layer. Containment of the PDR 
shell by the ICM prevents further penetration of FUV photons into the cold clump and there is no 
additional heating. For these clumps, photoevaporation is relatively unimportant and they retain large 
fractions of their initial masses. 

Clumps that are heated with t P uv > A evolve quasi-statically, steadily adjusting themselves to 
the pressure of the heated PDR shell. The evolution of these clumps can be determined from simple 
steady-state equilibrium considerations. In their final configuration, the clumps may retain a fraction 
of their initial cold gas, surrounded by a warm PDR envelope, in pressure equilibrium with the ICM. 
We consider two extremes in density profiles of the clump, a constant density clump and a truncated 
isothermal sphere ( n oc 1/r 2 ) in pressure equilibrium with the ICM. 

First, we solve for the evolution of a clump of constant density. From conservation of mass, 

n c0 r c0 3 = nj r f 3 + n PDR (rp DR - rj) (21) 

where n f and n PDR are the number densities in the final cold clump core and the PDR, and r c0 ,r ; , 
and r PDR are the radii of the initial dump, remnant clump, and the PDR envelope respectively. As 
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the final configuration is a remnant clump in pressure equilibrium with the PDR and the ICM, the 
final density of the core n/ = n c o, and 

n c oC c 2 +an c oC c 2 -\-f3n c0 c c 2 = UperCpbr +cmpDRC 2 +fi n caC 2 { n PDRl n c<sV (= Pressure in the ICM) (22) 

Further, we know that the column density through the PDR envelope is N 0 , as the remnant clump is 
shielded from heating. Therefore, 

flpDR^PDR ~ rj) = N 0 ( 23 ) 

Using equations (21) and (23) to eliminate updr and r pdr, equation (22) can be written as 

^ = (1 + a + 0)£-# 1_7 -<* ( 24 ) 


where 
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(25) 


and r)f = ricorj/No. For complete photodissociation rjj — 0, or £ — 7/ c o . For the standard case a 1, 
j3 = 1 and 7 = 4/3, and for 77 c0 > 1 , 12 > 1 , equation (24) can be used to define a critical column density 
for complete photoevaporation. 

r] cr ~ 0.48i/ /3 (26) 

All clumps with an initial column density lower than rj cr are thus eventually completely photodissociated. 
Alternately, a clump with an initial column density r] cr has to be exposed to an FUV field {v) greater 
than or equal to that given by equation (26) to be completely photodissociated. 


Figure 7 shows a plot of the final column density of a clump, r/ ; , as it is exposed to various FUV 
fields, measured by u, for a = 1, 0 = 1 and 7 = 4/3. The figure shows trajectories followed by the 
clump column density as the FUV field is increased, as might happen when the clump slowly emerges 
into the PDR region from the shielded interior of the molecular cloud. Note that the abscissa is i y 
and not the initial column density r/ c o- The outer column of cold clump gas Nq that is initially heated 
is thereby discounted, and as soon as the FUV field is “turned-on at time t — 0, t/j — t ] c 0 — 1. As 
the FUV field (u) increases, the dump begins to lose mass due to heating of the outer shell, and its 
column density begins to decrease with increasing v. As v approaches the critical value for complete 
photodissociation for that particular r/ c 0, the column density rapidly decreases with u and the clump 
photoevaporat.es completely. If the final strength of the FUV field is less than that required for 
complete photoevaporation, the clump column density remains fixed at a point on its trajectory, as 
determined by the value of the parameter v. Therefore, given a clump’s initial column density and the 
strength of the local FUV field, the final column density of the clump can be determined. Equivalently, 
the observed r} f and v can be used to determine the initial column r? c0 of the clump. The total mass 
loss of the clump during its evolution can also be evaluated from equation (24). Figure 7 also shows the 
contours for the mass loss being 5%, 50% and 95% of the initial mass, for different values of ?// and u. 


The evolution of a clump with an isothermal density profile can be determined analogously. The 
initial density distribution in the clump is given by n(r) — n cs or^ 0 /r 2 , where n cs 0 is the density at the 
clump surface, and the initial mass of the clump is 47r?rif/7r Ci o r co- ass ” rae that the density in the 
heated PDR shell in the final c onfiguration of the clump is constant. Utilizing equations( 21-23) for 
this case, we obtain an equation very analogous to equation (24) where now 




1/2 


( 27 ) 
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VcsO = n cs0 Tco/No and V f = n cs0 r f /N 0 . An effective average column density through the clump can be 
defined by the ratio of the mass to area of the clump, 7/ c0 = m c0 /(nrl 0 m H N 0 ) = iri cs0 • The evolution 
of a clump with an isothermal density profile is thus qualitatively similar to that of a constant density 
clump, but with a different effective column density for the same mass. These solutions apply to 
clumps that are massive enough for gravity to be important in determining their density structure, and 
therefore have appreciable central density concentrations. 


6. Numerical code and results 

In the above analysis, simplifying approximations were made, which are now relaxed in a numerical 
code to obtain complete and more exact solutions (see Appendix D for details). One of the major 
differences between the numerical code arid most of the analytic equations is the inclusion of gravity in 
the numerical code. 


6.1. Impulsively heated clump in a vacuum 

We first discuss the results of a numerical simulation of a clump generated by turbulence and 
heated impulsively. This case is hereafter referred to as Case V (for “vacuum”). A vacuum boundary 
condition is used so that the clump is free to expand, even in the absence of any heating. Initially, 
the thermal, turbulent and magnetic pressures in the clump are all assumed to be equal (a — — 1), 

and the magnetic pressure scales with the density to the power 7 = 4/3. The FUV field is turned on 
instantly (t FU v= 0)- The clump is assumed to be a constant density sphere, with a column density given 
by Vco = 10, and in a FUV field characterized by v - 10 in the outer layers. The value of t] cQ corresponds 
to that typical of clumps in clouds, for example, m c = 0.8 M e ,n c = 2 x 10 5 cm _3 ,and r c = 0.03pc. 
Clumps exposed to FUV fields with G 0 ~ 10 5 are heated to T PDR ~ 1000K on their surfaces; if their 
shielded centres are T c ~ 10 K, then the ratio of sound speeds u ~ 10. These initial conditions imply 
the clump lies in the shock compression regime of Figure 2, and we expect from our analytic solutions 
that it should photoevaporate in less than one sound crossing timescale t c (from eq. [11]) 

Figure 8 shows the density, velocity, pressure and temperature profiles in the gas as a function 
of radius at three different instants of time, t = 0.05t c ,0.1f c and 0.15/ c . As the outer layers heat 
up and the increased thermal pressure causes them to expand, a shock propagates into the clump, 
compressing the cold inner clump gas. The shock strengthens as it progresses to the centre of the cloud, 
as evidenced by the increasing density and velocity at the position of the shock front. The central 
cold gas gets compressed to very high densities, mainly due to the convergence of the radially-moving 
shocked gas. In Figure 8 we only show the evolution prior to the shock reaching the centre of the cloud. 
After the shock reaches the centre, the compressed gas rebounds and the cold inner clump undergoes 
radial oscillations as it settles into pressure equilibrium with the warm, expanding PDR outflow in the 
outer layers. The mean density of the compressed clump gas increases with time, as expected from our 
analytical solutions (eq. [B7]). As the clump loses mass, the FUV penetrates deeper and the clump is 
completely photoevaporated in 0.7 t c . Our analytical solution (eq. [11]) predicts a photoevaporation 
timescale of 0.5 t c , which agrees very well with the numerical solution. The photoevaporation timescale 
derived for a clump with this column density, r) c0 = 10, and density, radius and sound speed given by 
n c = 2 x 10 5 cm -3 , r c = 0.03pc,and c c = 0.3kms _1 respectively, is 10 4 years. 
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6.2. Gravitational collapse 

The compression of a clump by strong shocks due to FUV heating can raise the central densities 
by large factors. This, however, did not render the clump gravitationally unstable to collapse in the 
above case because 7 was chosen to be 4/3, which raises the magnetic pressure in dense gas sufficiently 
to prevent gravitational collapse (Chandrasekhar 1961). The scaling index of the magnetic pressure 
with density in clumps depends on the orientation of the magnetic field in the clump, and on whether 
the field is frozen into the neutral gas. For a frozen, uniform, unidirectional magnetic field, P B oc n 4/3 . 
The magnetic field configuration in clumps is considerably more complicated (Ward-Thompson et, al. 
2000), and ambipolar diffusion may allow the neutral gas to slip past magnetic flux lines; making it 
difficult to define 7 unambiguously. 

We investigate the possibility of triggering gravitational collapse in clumps through a simulation 
where we adopt 7=1, and all other parameters identical to those in Case V. Such a clump is expected 
to undergo gravitational collapse based on our earlier analysis (eq. [16] where ?/c0 — 10? an d Figure 
5). The evolution of the clump is initially similar to that in Case V. The outer surface expands, and 
a shock is driven into the clump centre. As the shock reaches the centre, the density of the cold 
clump gas increases and the clump radius decreases, until the central region becomes gravitationally 
unstable. A collapsing “core” is thus formed. This can be seen in Figure 9, which show's the mass of 
the clump as a function of time. The mass of the collapsing core is also shown in the figure, along w ith 
the results of Case V. The core mass is seen to increase above that contained in the clump just before 
collapse as the core accretes mass from the w'arm PDR gas. A significant fraction (60%) of the initial 
clump mass is driven gravitationally unstable to collapse, on very short timescales. In the absence 
of sufficient magnetic field support, exposure to strong FUV fields can thus trigger star formation in 
clumps. We also conducted a simulation wuth a weaker FUV field (v — 5), and found that the central 
regions collapsed in that case too, but the collapsing core contained only about 10% of the initial c lump 
mass. Although these triggered collapse solutions are interesting and instructive, we emphasize that 
for initially stable clumps with /3 ~ 1, they hold only for 7 < 4/3, w'hich requires ambipolar diffusion. 
Since ambipolar diffusion timescales are longer than the shock compression timescale, it is likely that 
realistic clumps require 7 > 4/3 and that triggered collapse does not often occur. 


6.3. FUV heating of a pressure-confined clump: tpuv t c 

Pressure confined clumps are modelled as Bonnor-Ebert spheres which are in hydrostatic 
equilibrium. We construct a Bonnor-Ebert sphere by setting up a density distribution for a clump 
wuth thermal, magnetic and turbulent pressure support and in hydrostatic equilibrium. The density 
in the clump decreases as r -2 in the outer regions (similar to an isothermal sphere distribution), and 
flattens out in the central regions. The ratio of central to surface density x, determines the stability 
of the sphere against gravitational collapse, with a value exceeding 14.4 being “critically unstable for 
an isothermal, non-magnetic: g<is. In these runs wuth a confining ICM, the outer 200 zones are used 
to represent the ICM. The ICM has the same pressure as the surface of the clump and is a constant 
density medium wuth a density (and temperature) contrast of 1000. As discussed in §4, for clumps 
confined by an ICM, a slow turn-on of the FUV field may be more appropriate. The FUV field in the 
numerical run is turned on in a time t FUV = 5 t c , and w'e increase the temperature on this timescale so 
that the heated outer shell reaches its maximum temperature at bt, c . 
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Figure 10 shows the results for a case with varying x and where 77 ^ = 10, v - 5,7 = 4/3, and 
a = P=l. For Bonnor-Ebert spheres with x > 6, the outer region with an isothermal density profile 
contains more than half the mass, and the clump evolution resembles that of an isothermal sphere. In 
the figure, the mass of the clump is shown with time along with the analytical results for an isothermal 
sphere and a constant density sphere with the same initial column. The numerical results for three 
values of the central to surface density ratios, x = 1,3 and 6 are indicated. For central concentrations 
greater than that represented by x = 6, the analytical results for an isothermal pressure-confined 
sphere match the numerical results, and the agreement increases as x increases and the Bonnor-Ebert 
clump configuration approaches a 1/r 2 density profile. For x = 1, the clump is a constant density 
sphere, and the corresponding analytical calculations apply. A slight discrepancy can be noted in the 
figure between the two results, and this is due to the differences in the more realistic numerical model, 
i.e. an exponential drop in the PDR shell temperature with column and the inclusion of gravity . For 
intermediate central density concentrations, as represented by the case with x = 3, the mass loss rate 
is in between the two analytical extremes, as expected. 


7. Acceleration of clumps due to rocket effect 


Spherically asymmetric m.iss flows from photoevaporating clumps can cause a rocket effect on 
the clump, as first noted by Oort & Spitzer (1955) in their theory on cloud acceleration by ionizing 
radiation from OB stars. Clumps exposed to FUV radiation in PDRs see a stronger incident flux of 
photons from the direction of the OB stars, though a significant flux due to backscattered radiation 
from dust on the opposite surface tends to reduce this asymmetry (see §5.2). Ihe front and back 
hemispheres of the clump (with respect to the OB star) are thus heated to different temperatures 
(different v) and the escaping material has a higher flow velocity on the star-facing surface of the 
clump. The net thrust on the clump accelerates it away from the star, causing the rocket effect. With 
the help of our analytical expressions for the evolution of an impulsively-heated, photoevaporating 
clump in a PDR, we now evaluate a simple criterion which determines the clump size at which the 
rocket effect becomes significant. 

A simple 1 -dimensional formulation for the acceleration due to the rocket effect is given by 


dVR 

dt 


1 dm c 
m c dt 


(v ff - v fb ) 


(28) 


where v ff and v fb are the flow velocities at the front and back surfaces respectively and are assumed 
equal to the PDR sound speed Cpdr at those surfaces. We first compute the acceleration due to the 
flow from each hemisphere separately, take the difference to estimate the net acceleration, and then 
calculate the net velocity v R attained by the clump. The mass loss from a hemisphere is given by 
2nr c (t) 2 pbCp D R, where p b is the density at the base of the flow. Taking only the component of the 
momentum in the flow along the direction to the star, we obtain 


f~rO = -^-^-nr c (t) 2 p b cp D R (29) 

V dt y 1/2 m c 

The net acceleration on the clump is therefore given by 

~j ~ \(Pb c2 pDR)f ~ ( PbCp DR )b } 
dt m c 


( 30 ) 



- 20 - 


where the suffixes / and b denote the quantities for the front and back surfaces of the clump, 
and the clump radius 77 c is determined by the dominant flux on the front surface. As the initial 
expansion/compression phases in the evolution of a photoevaporating clump are typically short (see 
§5), we ignore this initial phase of evolution and assume the clump is evolving such that its column 
density is at its critical value 77, T ,;* ~ 4^ 2 /3, for the standard values of o, /? and 7 (see eq. [6]). Along 
with the relation for the column through the flow nb(t)r c (t) = 2 Nq for each surface and the mass of the 
clump m c = 4nmHTif c r 3 j c /3, (see appendix, eq. [A2] ) , we thus obtain for the net acceleration of the 

clump, 


dvi 


9 cl 


dt 8 r. 
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PDRJb 


( c pdr )/ J 


(31) 


The clump will remain in the PDR of the GMC and be subject to photoevaporation if the rocket 
velocity v R attained in a time tpo R = X pdr/ v if is l ess than vj F , or equivalently if 


9 cl 
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8 v 2 
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( C PDJ ?)/ 




PDR 


(32) 


There is thus a critical size r R of clumps, where for r c < r R they get rocketed back into the shielded 
molecular cloud, and for r c > r R remain in the PDR and get photoevaporated. For example, for front 
and back surface flow velocities of 3 kms -1 and 2 kms -1 respectively, equation (32) gives a critical 
clump size 

r R = 0.225A t -pdr (33) 


Clumps much smaller than this will get rocketed back into the molecular cloud as soon as they enter the 
PDR, partially evaporating in the process. We emphasize again that the rocket effect is significant only 
for turbulent clumps which undergo impulsive heating and photoevaporative flows. Pressure-confined 
clumps are gradually heated on their exteriors and do not go through phases of evolution with rapid 
mass outflows, and hence experience an insignificant rocket effect. 


The mass lost by a clump before it gets accelerated out of the PDR. can be estimated by integrating 
equation (28) with time, where only the components of velocity and mass flow along the direction 
towards the source are to be considered. Thus the mass of the clump is approximately given by 


m c (t) = m c oexp(—2vjp/ A cpd R ) (34) 

where A c P dr is the difference tn the flow velocities on the two sides of the clump. For typical values 
of A c PDR = 1 km s -1 , and v, F = 0.5 km s _1 , in the time it takes for v R to reach v IF , the clump loses 
about 63% of its initial mass. Therefore clumps which are rocketed back into the cold molecular cloud 
interior also lose significant fractions of their mass before they attain high enough velocities to keep up 
with the advancing ionization front. 


8. Discussion 

We begin our discussion by summarizing our results for both the case where clumps are produced 
by turbulence in regions with insignificant interclump pressure and the case where clumps are initially 
stable structures confined by interclump pressure. 

Consider first the scenario where clumps are continually generated by turbulence in the molecular 
clouds and the PDR. Initially, i,he evolution of the clump is determined by its column and the strength 
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of the FUV flux. If A > 1, or r; c0 > 4iz 2 /3 for a = 1, 0 = 1 and 7 = 4/3, there is a very brief period 
of shock compression but the shock quickly dissipates before making it to the clump centre and the 
clump radius shrinks only slightly. Subsequently, the clump expands until Aw 1, at which point 
the photoevaporative flow halts the expansion and proceeds to shrink and compress the clump gas, 
maintaining A = 1, or a column n c r c = 4v 2 N 0 /3. If A < 1, then the shock propagates to the centre of 
the clump, compressing it significantly, followed by an evolution at A ~ 1 identical to that above. 

During the final A = 1 evolution of the shrinking, compressing, evaporating clump, the rocket 
effect may play a significant roie depending on the physical size of the clump at this time. If r c ^ r R , 
the clumps wiil be accelerated into the molecular cloud, losing a significant fraction of their mass in the 
process. If r R <r c < X PDR , the clumps will tend to evaporate in the PDR region. However, because 
in the turbulent model clumps are forming in the PDR, including the region near the IF, some of 
these clumps will survive to enter an advancing H II region. If r c > X P dr and a blister H II region is 
advancing, then most of these large clumps will survive to be engulfed by the H II region. 

The lifetimes of turbulent dumps are of order several sound crossing times of the initial size of the 
clump (Eq. [1 1] and Eq.[20]). As discussed previously, this is somewhat longer than clump lifetimes in 
the absence of an FUV field. Since turbulent clumps are being continuously formed, the steady-state 
abundance of clumps could actually increase in the presence of an FUV field. 

Based on these results, the following predictions can be made for the turbulent scenario of 
photoevaporating clumps in PDRs. 

1. There will be an enhanced population of clumps with columns n c r c w 4v 2 N 0 /3. Note that this 
implies that in comparing this population from one PDR to another, the PDR with the lower u 
(e.g., due to lower G 0 ) will have clumps with smaller column densities. 

2. There will be a population of small clumps ( r c < r R ) moving with velocity v ~ vj P into the 
molecular cloud at A v > l. 

3. Compared with clumps in FUV-shielded regions, the clumps in turbulent PDRs will be smaller, 
denser and potentially more numerous. 

4. Similarly, the clumps entering the H II regions will be smaller and denser than the molecular 
cloud clumps, and will have columns peaked at 4/AV 0 /3. 

Consider next the scenario where clumps are pressure-confined structures which are bounded 
by a low-density ICM, and the FUV turn-on timescale is of order the crossing time X PDR /v, F for 
clumps through the PDR. For small clumps with r c < X PDR c c /vj F , the FUV turn-on time t F uv > t c , 
which means that the FUV quasi-statically heats and expands an outer PDR surface on the clump, 
continually maintaining a pressure Pjcm in the PDR region and in the cold clump region. Again, the 
evolution of the clump is determined by its column and the strength of the FUV flux. If ry c0 < 0.48A 3 
(for a = l,/3 = 1, and 7 > 1), the clump is completely heated and photodissoeiated as it enters the 
PDR from the molecular cloud. If r) c0 > 0.48;/ /3 , the clump shrinks to a smaller cold clump surrounded 
by a PDR shell. Since the evolution is quasi-static, there is insignificant rocket effect on these clumps. 
In this scenario, small cold clumps, protected by PDR shells, can survive passage through the PDR 
and into the H II region. However, the minimum initial column for such a clump is rj c 0 ~ 0.48A 3 (for 
a = 1, /3 = 1, and 7 > 1). 
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For large clumps, r„ > X P1 , R c c /v IF , l FU v < <„ which means that shocks compress the clump 
and strong photoevaporative outflows are initiated as the clump enters the PDR Clumps \u • 

1 < „ o < 4i/ 2 /3, are compressed by shocks which propagate to the centre, but they finally re ax 
to an analogous situation as above, i.e., a thick PDR shell surrounding a small mass cold clump if 
Vc0 > 0.48i/ 4 / 3 and totally photodissociated otherwise. Clumps with 77 c0 > 4i/ /3 have a brief period of 
shock compression and PDR expansion, but rapidly relax to a configuration with a thin PDR shell and 
most of the initial mass still in the cold shielded core. All the large clumps experience a rocket effect 
However, as shown earlier, large clumps with r c > r R do not experience enough rocket acceleration to 
prevent them from entering the H II region. Therefore, in the ICM confined scenario, a variety of clump 
sizes enter the H II region. However, the small clumps entering the H II region were initially clumps o 
larger size and mass in the molecular cloud, and they enter the H II region with large protective PDR 

shells. 

One can compare the abundance of clumps of various sizes in PDRs with the abundance of clumps 
in the molecular cloud, under the assumption that ICM-confined clumps are longdived (the formation 
time of clumps exceeds the crossing time X PDR /v [F ). Small cold clumps are converted to PDRs but 
larger cold clumps are converted to small cold clumps with PDR shells and this partially replaces the 
small clump population. Thus, there will be a significant drop in the cold clump population with 
n c0 < 4^/3. However, the population ^ > 4^/3 is little affected. Therefore, a variety of clump sizes 
enter the H II region, but there is a suppression of clumps with columns i** < 4^ 2 /3 compared to the 
distribution in the shielded molecular cloud. 

Based on these results, the following predictions can be made for clumps initially confined by 
interclump pressure. 


1. Compared with small cold clumps (r/ c0 < 4i/ 2 /3) in similar pressure regions of the molecular 
cloud, clumps of the same column density in PDRs will have a smaller relative population. 
Correspondingly, small and large clumps enter the advancing H II regions, but the relative 
number of small clumps is suppressed. 

2. Larger clumps are relatively unaffected, but do experience an initial transient period of shock 
compression and photoevaporative flow. 

3. Small chimps ( Vc0 < 4iz 2 /3) in the FUV-illuminated region will have small masses relative to 
their PDR shells. 


We have shown that small clumps tend to be destroyed before they are overtaken by the ionization 
front of the advancing H II region. This would affect the number, size distribution, and structure 
of pressure-confined clumps that enter an H II region, and therefore affect the evolution of the H II 
region (see Bertoldi 1989 and Bertoldi k McKee 1990 for a discussion of the effect of photoevaporating 
clumps on the evolution of an H II region). Clumps that survive to enter the H II region will thus 
have undergone significant evolution from their initial states. These clumps could be considerably 
denser if they are still vigorously photoevaporating as they enter the H II region, where they continue 
to photoevaporate (Bertoldi 1989, Bertoldi & McKee 1990). We will discuss the propagation of clumps 
into H II regions in a separate paper, with application to the Eagle Nebula. 

We conclude our discussion by using our results to critique the previous inferences of dense (~ 10 
enr 3 ) FUV heated clumps in PDRs and to propose several observational consequences of our results. 
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Our photoevaporation models provide arguments against the interpretation that dense (n ~ 10 
cm -3 ) PDR surfaces of clumps give rise to the observed high excitation CO lines. The fundamental ^ 
problem is the propagation of the FUV photons to regions of such high density. If densities of ~ 10 
cm" 3 are in fact heated to T » 100K, then the thermal pressures in these regions are much larger 
than the pressures in the interdump gas. The only way this is possible is if: (i) the clumps are 
gravitationally bound or (ii) the warm surface regions are photoevaporating, as discussed in this paper. 
In either case, the density in the clump will smoothly fall off with radius, n(r) oc r where b ~ 2 - 3, 
until the transition to the ICM. However, FUV photons which heat the CO are only able to penetrate 
through a column No ~ 2 x 10 21 cm -2 , where 


N o 



n(r)r 

b- 1 V ' 


(35) 


This equation reveals the size r of the FUV heated region if a density n(r) is required in this region. It 
implies that for the FUV to heat gas at densities > 10 6 cm 3 , the clumps must be very small with size 
scales of r ~ Vo/10 6 cm -3 ~ 2 x 10 15 cm at this density. Gravitationally bound clumps have densities 
n ~ 10 6 - 7 cm" 3 at sizes » 10 15 cm (e.g., Shu 1977). In such clumps, the FUV would not penetrate to 
such high densities, but would be absorbed in the lower density regions further out. Such small scales 
as 10 15 cm also argue against unbound clumps explaining the high temperature CO emission. We have 
shown in this paper that such a small clump will either be rocketed out of the PDR or photoevaporate 
in roughly a sound crossing time, or ~ 300 years, a time so short as to rule out this possibility . 

The size scale of 100 AU is suggestive of disks around young stars rather than clumps, so that 
one might appeal to photoevaporating disks like the “proplyds” in Orion (e.g., Johnstone et al. 1998). 
Here, the enormous mass reservoir of the disk greatly lengthens the photoevaporative timescale. 
However, the area filling factor of these photoevaporating disks required to match the CO high-J line 
intensities in, for example, the Orion Bar region is of order ~ 0.1 (Burton et al. 1990). Such a high 
area filling factor requires a volume density of > 3 x 10° stars pc 3 , more than an order of magnitude 
higher than that observed in the densest region of the stellar cluster near the Trapezium (McCaughrean 
et al. 1994). Photoevaporating proplyds are thereby also ruled out. 

Therefore, it seems quite unlikely that FUV is heating dense, n ~ 10 6- ' cin" 3 , PDR surfaces of 
clumps or that photoevaporating proplyds can explain the mid-J CO emission. We suspect that the CO 
emission arises in n ~ 10 5 cm- 3 surfaces of large (> 10 16 cm) clumps, and that the chemistry, heating 
and/or dynamics needs to be modified in PDR models in order to better match the observations. 

There are several observational consequences implicit in the scenario of a strongly photoevaporating 
(i.e., tpuv/tc < 1) clump. First, the molecular material from the clump core is effectively advected 
out into the PDR region, as has been discussed by Bertoldi k Draine (1996), Storzer k Hollenbach 
(1998), Hollenbach k Tielens (1999). This can have strong effects on the chemistry, primarily because 
H 2 can exist in surface regions where it would have been atomic in a stationary, steady-state case. 

In turn, the higher abundance of H 2 can chemically enhance the abundance of other minor species, 
such as H^ (Bertoldi k Draine 1996) and CO and can modify the heating and cooling processes and 
therefore the PDR temperature. Initial results (Storzer k Hollenbach 1998 and Storzer 2000, private 
communication) suggest that advection will enhance the intensities of the mid-J CO lines, by creating 
more warm CO near the surface. Finally, the flow will broaden line widths to values A vfwhm ~ c pdr- 
However, it should be noted that many common molecules, such as CO, HCO + , CS, HCN, CN still 
exist primarily at column depths N > N 0 from the surface, where the advection velocities are low. 
Broader linewidths should be seen in C ll(158//m),0 1(63, mi) and possibly the H 2 (2pm) lines. 
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9. Applications to observed PDRs in well-studied star forming regions 

We compare our results with available observational data on clump characteristics from some 
well-studied star-forming regions, the Orion Bar, M17SW, NGC 2023 and the PDR surrounding the 
Rosette Nebula. Clump column densities and hence 77 can be estimated from measured or inferred 
densities and sizes of clumps in these PDRs. Densities are usually inferred indirectly from comparing 
observed line intensities and line ratios with available PDR chemical models in most of these cases. 
Clump sizes are ill-determined, as they are usually too small to be resolved directly and only upper 
limits to the sizes are available with the PDRs in Rosette and NGC 2023 being possible exceptions. 
The FUV field strengths in those PDRs can be estimated more accurately, as the luminosities of 
the illuminating O/B stars is reasonably well-determined. We use these field strengths and the dust 
temperature models of Hollenbach, Takahashi & Tielens (1991) to estimate the temperature of dust in 
the shielded clump interior. As the interior gas temperature is closely coupled to the dust temperature, 
we set the two equal and thus determine the sound speed c c in the cold clump gas. The temperature 
to which the FUV field heats the surface clump gas is determined from the PDR models of Kaufman 
et al. (1999), and thus c, PDR is obtained, allowing us to estimate the value of the parameter u in each 
region. Below we first determine 77 and v for each observed region, and then use these parameters, 
compared to Figures 2 and 7, to infer their nature and evolution. 


Orion Bar The Orion Bar h;is been observed in many atomic and molecular line transitions in the 
sub-millimeter and infrared wavelength regions (e.g. van der Werf et al. 1996, Tauber et al. 1994, 
Young Owl et al. 2000, Marconi et al. 1998). The molecular line emission peaks at about 0.03 pc 
from the IF, and the gas between this dense ridge and the IF is mainly neutral atomic hydrogen. The 
observed line intensities and line ratios appear to be best explained by a two-component model for the 
PDR gas, with small (< 0.02 pc) dense (n c ~ 10 6 " 7 cm" 3 ) clumps, required to match the line emission 
and an interclump medium tijcm ~ 10 4 " 3 cm -3 which causes the chemical stratification of the edge-on 
PDR. The local FUV field is estimated as G 0 ~ 4 x 10 4 . Such an FUV field incident on clumps with 
these inferred densities would heat their surfaces to temperatures ~ 2000 - 5000 K (Kaufman et al. 
1999). The gas temperature in the interior of clumps exposed to this FUV flux is expected to be ~ 50 
K (Hollenbach et al. 1991). We thus obtain 77 ~ 15 — 150 and v ~ 8 — 10. 


M17 SW Multilevel molecular line observations in CS (Wang et al. 1993, Evans et al. 1987) and in 
several fine structure lines (C II, Si II, O I, Meixner et al. 1992) indicate that the M17 SW cloud core 
consists of numerous high density clumps with densities n c ~ 1 — 5 x 10 5 cm 3 . Maps of CO + emission 
compared with theoretical models suggest that the emission is produced in the warm surface layers of 
PDRs in dense clumps (Storzer et al. 1995) of sizes r c ~ 0.1 pc. The local FUV held near the H II 
region/molecular cloud interface is Co ~ 8 x 10 4 , which implies PDR temperatures of about 2000K 
on the clump surfaces and cold clump gas temperatures ~ 60K, therefore, u ~ 4 — 9. The column 
densities of clumps are typically 77 = 7 — 35, as estimated from typical densities and sizes (Meixner et 
al. 1992, Storzer et al. 1995). 


NGC 2023 The well-studied PDR in the reflection nebula NGC 2023, at a distance of 475 pc, 
is illuminated by a B star and is well-studied (e.g., Wyrowski et al. 2000, Steiman-Cameron et al. 
1997, Howe et al. 1991). Observations of atomic fine structure lines indicate the presence of clumpy 
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gas in the PDR with densities ~ 10 5 cnT 3 , and an ICM density of ~ 10 3 cm" 3 (Steiraan- Cameron 
et al. 1997). The PDR is situated about 0.1-0. 2 pc away from the exciting star, and has a width 
Xpdr = 0.04 pc. The clumps have sizes ~ 0.05 pc (Wyrowski et al. 2000), and therefore r) c0 ~ 5. 
The FUV field strength has been estimated at G 0 = 10 4 , and from the PDR models of Kaufman et 
al. (1999), c PDR = 3kms -3 . With a clump gas temperature of ~ 20 K for this value of G 0 , we obtain 
i/ = 5. It should be noted that the levels of turbulent energy in the NGC 2023 PDR as estimated 
from observed linewidths are low (Wyrowski et al. 2000), which suggests that clumps are not rapidlj 
re-formed by turbulence in the PDR. 


Rosette PDR The Rosette Nebula surrounding the young open cluster NGC 2244, is an expanding 
H II region at a distance of 1.6 kpc. The H II region/molecular cloud interface shows a distinct ridge 
of emission, with substructure on small scales down to 0.5 pc (Schneider et al. 1998). The PDR. region 
is very extended, and the FUV photons penetrate deep into the cloud with an interclump medium 
density n ICM ~ 10 cm" 3 (Blitz 1991, Williams et al. 1995). The ICM near the edge of the HII region 
may even by partially ionized gas as suggested by the pervasive Ha emission observed (Block et al. 
1992). This region has a low FUV field, with Go ~ 200 about 15 pc away from the OB star cluster, 
and dropping to lower values of 10-50, about 30 pc away from the cluster (Schneider et al. 1998). 
Clump densities required to match observed CII line intensities are ~ 10 4 5 cm 3 . Cold clump gas 
temperatures are calculated to range from 10 — 15 K. The Rosette PDR thus has r/ c0 ~ 5 — 35 and 
i/ ~ 5, while further away from the H II region, u ~ 2 — 3. 


Using the above determined values of the parameters rj and v for clumps in the above PDRs, we 
can locate them on the parameter plots of Figures 2 and 7. Figure 11 shows an overlay of our model 
results with inferred data for both the turbulent and pressure-confined clump models. The data are 
consistent with both clump models and, within the limits of the uncertainties in arriving at “observed” 
rj and v, agree reasonably well with our predictions. Turbulent clumps are expected to have columns 
fairly close to the critical value i] C rit which is ~ 4^ 2 /3. The observed data lie near the r] = rj cr n line, 
indicating that the observed clumps may be turbulent in origin and undergoing photoevaporation in 
the FUV field. The exception is NGC 2023. The clumps in this PDR. are observed to have negligible 
turbulence, and are thus not expected to have a column equal to Tj cr n . The data are also consistent 
with clumps being pressure-confined. In such a situation, we would not expect to see clumps beyond 
the 95% mass loss line, as seen in the figure. However, there does appear to be a weak trend of 
increasing t/ with u. which if real is not easily explained in this model. The clumps in NGC 2023 are 
probably pressure-confined, and seem to have lost almost 75% of their original mass through FUV 
heating. High resolution observations of clumps in PDRs, with more accurate determinations of the 
sizes and densities of clumps would help in locating them more precisely on the rj — v plots, and thus 
be able to distinguish between the pressure-confined and turbulent clump models. 


10. Conclusion 

We have studied the effects of FUV radiation from young OB stars on the evolution of dense 
structures or clumps in photodissociation regions. We find that the ambient FUV field penetrates 
through the surface of a dense clump and heats this surface layer to high temperatures, causing mass 
loss and thereby inducing pholoevaporation of the clump. Through analytic approximations and 
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numerical hydrodynamical calculations, we determined the evolution and lifetimes of clumps, subject 
to various physical conditions and initial parameters. 

The evolution of a clump is mainly determined by three parameters, the ratio of the initial column 
density to the column penetrated by the FUV flux, t? c0 ; the strength of the FUV field, which we denote 
by the ratio of the sound speeds in the heated gas and the cold clump material, v\ and the timescale 
for the turn-on of the FUV field, relative to the sound crossing time through the clump, t FUV /t c . We 
consider the evolution of turbulent, impulsively-heated clumps and also pressure-confined clumps in a 
PDR. Impulsively heated clumps with a photoevaporation parameter A « 3r/ c0 /(4i/ 2 ) < 1 are initially 
compressed by a shock driven by the external high-pressure FUV -heated layer to high densities, and 
later lose mass as the FUV gradually penetrates through the entire clump. Clumps with A > 1 go 
through an initial expansion phase after the rapid decay of the shock wave. Both clumps with A < 1 
and A > 1 evolve toward the A = 1 condition. The final evolution toward complete photodissociation 
and heating occurs at constant column through the clump, the density increasing as the clump shrinks. 
Photoevaporation timescales in all these cases are typically a few t c , suggesting clump destruction 
times of ~ 10 4 ~ 5 years, under typical PDR conditions. Slowly heated clumps evolve quasi-statically, 
and if confined by an external I CM pressure, a fraction of the initial cold clump mass may be retained 
if ijco > 0.48iA Clumps with initial column densities r ) c 0 < 0.48// 2 are completely photoevaporated by 
the FUV flux. 

We predict that in the turbulent scenario, observed clumps should all have columns close to the 
critical value for the local FUV field, r] crit « 4^ 2 /3. Clumps lifetimes are prolonged by photoevaporation 
and we also expect a higher steady-state abundance of clumps in PDRs as compared to the shielded 
molecular cloud interior. For pressure-confined clumps, there will be a decrease in the number of 
smaller clumps in PDRs, and many clumps will lose substantial fractions of their mass t.o extended 
warm PDR shells around them We compared our results to observations of some well-studied PDRS, 
such as the Orion Bar, M17SW. NGC 2023 and the Rosette PDR, and find that the data are consistent 
with both clump models, but perhaps favour the turbulent clump interpretation. Clumps entering the 
H II region are expected to have thus undergone significant evolution as they pass through the PDR. 


We acknowledge support from the NASA Astrophysical Theory. Program under RTOP 344-04- 
10-02. U. Gorti was supposed by a National Research Council Associateship. We thank K. R. Bell, 
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A. Analytic model: Unconfined clumps with t F uy/t. c <C 1 

W T e first define several characteristic timescales during evolution. The initial sound crossing time 
through the clump is denoted by t c = t c q/c c . We can similarly define a sound crossing time through 
the warm PDR layer, as t p = 6«/c PD r. Finally, we define the timescale for the PDR layer to expand to 
a thickness comparable to r c0 , Go = Tco/c-pdr- 

W ; e relate the density n b at the base of the PDR shell (r = r c ) to N 0 for epochs when the PDR shell 
has expanded to a thickness greater than the clump radius ( t > Go)- The FUV '-heated shell of initial 
thickness <5o expands as an isothermal flow into vacuum and the characteristics of the flow regions are 
essentially similar to that of the Parker wind solution with no gravity (Parker 1958). For simplicity, 
it is assumed here that once the warm PDR region has expanded to <5 > r c or t > Go, the gas leaves 
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the clump surface at c. PDR and r,he velocity of the heated gas increases linearly with radius. Based on 
our hydrodynamical code results, this is a reasonable approximation to within about two clump radii, 
where most of the column density in the PDR gas lies. The clump mass loss rate is determined by 
the penetration depth of the FUV flux, and hence N 0 . For a power law density profile in the shell 
(n(r) oc r~ b ,b> 2), most of this column arises from the base of the shell, or, more precisely, from the 
region between the base (clump surface) r c and ~ 2 r c . From the equation of continuity and v oar, the 
steady flow density profile of the warm PDR shell is obtained as, 

n(r) oc u -1 r -2 oc r~ 3 , (Al) 

valid for r c < r < 2 r c once the shell has expanded to at least two clump radii (f > t e o)- Equation (Al) 
can be easily integrated over r to obtain the column through the warm PDR at any given instant of 
time. Let r c denote the instantaneous cold clump radius and r t the total radius of cold clump and 
warm PDR gas. At t, = 0, r c = r c (0) = r c0 - 8 0 . Using the fact that r t » r c for t > t e 0 , and as the FUV 
flux always penetrates through a column Nq. we have 

f Tt 1 

A 0 = / n(r , t)dr « -n b (t)r c (t) for t > t e0 (A2) 

J r c ~ 

We next write down the equations of momentum flux or pressure. The pressure in the shocked 
cold clump gas is the sum of thermal, turbulent and magnetic pressures and given by 

P c (t) = n c (t)c 2 c + an c (t)c 2 c + pn c0 c 2 r {n c {t)/n c0 y (A3) 

where n c (t) is the shocked clump density at time t , n c o is the initial clump density, a is the initial ratio 
of turbulent and thermal pressures, and /? is the initial ratio of magnetic to thermal pressure in the 
clump. The pressure at the base of the PDR flow, P b , is given by 

Pb(t) = n b (t){c 2 PDR + ac 2 c ) + pn c0 c 2 c + n b (t)v 2 lcmr (A4) 

where v f i ow is the outward flow at the base. We assume that vji ow = c PDR , and from equations (A2) 
and (A4), 

•Ffc(f) = 2 7 Tr( 2 cp£) R + ac . c ) + 3n c oc c 
rc{t) 

The boundary between the “shock compression” and “initial expansion” regions of Figure 2 
is determined by comparing equations (A3) and (A5). Values of the parameters t] c q and v on this 
boundary are such that the shock dies out just before making it to the centre. The pressure in the 
clump at this instant is about the same as the initial pressure and r c (f e o) — r c (0). In a time t = t e o, 
the outer edge of the PDR shell expands out to 2r c0 , and the initial PDR column begins to decrease 
appreciably due to spherical divergence, allowing further FUV flux penetration. 

If p b (t e 0 ) > P c {t = 0), the pressure in the PDR flow is sufficiently high for the shock wave 
to reach the centre of the clump, compressing the clump and raising its central pressure so that 
p b (t e 0 ) = P c (t e0 ). Such clumps lie in the “shock compression” of Figure 2. If the initial parameters are 
such that P b {t e 0 ) < P c (t, = 0), the shock wave rapidly dies out before reaching the centre and the clump 
then stops contracting and begins to expand at a speed « c c , maintaining pressure equilibrium. These 


n t 


2Nq V 
:0 r c (t)J 


(A5) 



two evolutionary sequences are separated by the (7^0, v) boundary obtained setting Pb(t e 0) — Pc(t — 0), 
or from equations (A3) and (A5). 


Vco ~ 1 \ VcO 1 / 


1 +« + /?. 


(A6) 


For values of 7/ c0 and v typical in PDRs such as Orion (r/ c0 > 1,/A > l,a = 0 = 1,7 = 4/3) this 
relation can be approximated as 77 c o ~ 4/A/3. A photoevaporation parameter A = 3(7 c o — l)/2(2/A -t- 1) 
can be defined such that the relation A = 1 demarcates the two regions. Clumps with A < 1 are shock 
compressed and those with A > 1 initially expand. 


B. Evolution of shock-compressed clumps 


As the flow expands in a shock compressed clump, the clump shrinks by losing mass at a rate 

^ = -4tt m H n„(t)r e (t) 2 cp DR (Bl) 

at 

Initially, a shock is driven into the clump, but we do not attempt to analytically model the propagation 
of the shock wave in detail. We assume that the shock travels inward at the PDR sound speed and the 
entire clump gets compressed in a time t s = r c (0)/cpDR , where r c (0) = r c0 — <5 0 . We first solve for the 
clump radius, mass and mass loss rate for t < t s , during the epoch when the shock propagates to the 
centre. 

We assume that the clump radius decreases with a constant speed cppp) so that 


r c (t) = r c ( 0 ) - v b t. 

The mass of the cold clump decreases with time, and at t s is given by 


m, 


.( t s ) = m c (0) + f ( dm c /dt)dt 

Jo 


(B2) 


(B3) 


where m c (0) = |7rm//n c0 r r 3 (0). From equations (A2), (Bl), (B2) and (B3), we have 

m c (t t ) = m c (0) — 8T:m H Nor c (0) 2 (l — v b j{2c P p,R )) (B4) 


We solve for v b and thereby calculate r c {t s ). To do this we need to calculate the clump compression 
at t s . The clump gets compressed by the shock to a new density n(t s ) which can be determined from 
the condition of pressure equilibrium at the clump surface. The momentum flux conservation equation 
across the “front” which marks the boundary between clump and FUV heated flow is then 

(1 + a)n c (t s )c 2 + fin c oc 2 (n c (t s ) /n c o ) 7 — 2 n b (t s )cp DR + anh(t s )c c + P n cO C c ( n b(ts) /ftco) 1 (Bf>) 

In moderate to strong FUV fields, the PDR sound speeds and hence the shock velocities (assumed 
equal to cpp>p) are high, and the resulting compression is also high. As the magnetic pressure Pp, 
scales with a higher power of the density than the thermal pressure Pp, (we use 7 > 4/3) and as 
the initial magnetic pressure is comparable to the thermal pressure, the thermal contribution to the 
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pressure in the compressed gas can be ignored. In the expanded flow, the density is low and here 
the thermal pressure and dynamical pressure dominate. Equation (B5) can be simplified to give the 
density of the compressed gas 


n c (t s ) ( (2;A + a) n b {t s ) 


^cO 


0 


n c o 


1/7 


(B6) 


From equations (A2) and (B6), 


n 


h(t s ) _ ( 2£p(2^ 2 + Q) y /7 r ^fJ-V7 


^cO 


The mass of the clump at time t s can also be expressed as 


TTl c {t s ) — ^ TTliiTl c (t s }r c (t s ) ■ 


(B7) 


(B8) 


From equations (B2), (B7) and (B8) 


4tt 

m c (t s ) = — m H n c o 


/2^(2^+Q)y /7 rc(Q) 3-i/7 ^ _ 


Vb 

CPDR 


3-1/7 


(B9) 


At time t„, the radius of the shock-compressed clump is r s - r c (t s ) - r c (0)(l v b /c PD fi) from 
equation (B2). The unknown parameter v b is finally determined from equating the two expressions 
for the clump mass, equations (B4) and (B9). The evolution of the clump radius is thus known as 
a function of time for 0 < t < t s . from equations (B2),(B7) and the solution for v b . These general 
equations are, however, fairly complex, and an approximate solution can be obtained for the case 
a = 1,0= 1, and 7 = 4/3 as, 

r, » ,r c (0) ( 1 - Hl+il'l (BIO) 


T]c 0 - 1 

where q = (A/3) 1 / 3 . 

We now solve for the evolution of the clump for t > t s . After being compressed by the shock, the 
clump continues to lose mass through photoevaporation, with a mass loss rate, 


dm c 

dt 


= -4irm, H n b {t)r c {t) 2 cp DR = - &nm H N 0 r c {t)c.p DR 


(Bll) 


The mass of the clump at any given time ( t > t s ) can also be written as 

m c t) = m c (0) 


2(2P + a )\' h (r c (t)' 3 -' h 


PiVcO - 1) 


^c(O). 


(B12) 


From equation (Bll) and differentiating equation (B12) with time, the radius of the clump as a 
function of time is given 


T-c{t) 


(=») 


2 1/7 (27-I) 6t/?7 e o / 0(^co ~ 1) \ 1/7 ( L - ^ 

(37 — 1) (wo — l) 2 \ 2 ( 2 iA -I- ct) ) \t c t] c qU 


(B13) 
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A photoevaporation timescale tpp for the clump can be defined as the time for the radius of the clump 
to shrink to zero, and from equation (B13) 


tpE — A 


r. \ 2 1/7 /3y - 1 ^ (t ? c0 - l) 2 ( 2(2zA + a) \ 1/7 + r/ c0 - 1 

VcOV 


T-c(O) 


27-1/ 6 vr ) c0 \ P{r] c o ~ 1 ) 


(B14) 


For a = l,/? = 1, and 7 = 4/3, the above equations determining the evolution of the clump can be 
simplified as 


and 


/.X II r s \ ^ IOtJcoZZ 9/4 / i (VcO 1) 

r<(<) = 1 U(0 ) ) "aSTiF* 


4/5 


r c (0), 




9/4 


tpE — tc 


m c (t) = m c { 0) 

3 ( 77 cq - l) 2 A _ 3(]_+9)y /y + (jjco ~ 1) \ 
lOgT/cQjz V %o - 1 / ^co*' / 


^c(O) / 

5/9 


(B15) 

(B16) 

(B17) 


C. Evolution of clumps that undergo an initial expansion 

Clumps with large initial column densities and in low FUV fields quickly evolve to the point where 
their internal pressures are greater than that of the (expanded) thin heated surface layer. They expand 
out into vacuum (as they would even in the absence of an external FUV field) at their sound speed c c 
for a time t e , until the pressure drops to that in the heated outer layer. At t = t e , there is pressure 
equilibrium and equation (B5) again holds. The clump expands to a new radius r c (t. e ) = r c (0) 4- c c t e 
and during expansion loses mass at a rate given by equation (Bl). The mass of the clump at t = t e is 
obtained as earlier, 

m c (t e ) = m c (0) - Snm H NoCpDR(r c (0)t e + c c t 2 e / 2). (Cl) 


Also, 


4 7T 

m c (t e ) = — m//u c (te)r c 3 (t e ) 3 - 


(C2) 


Because there is significant expansion, n c (t e ) <C n c 0, we assume that thermal pressure dominates in 
both the clump and the PDR flow. From equations (B5) and (C2) 


. . 87r 2 u 2 + a 2 

rn c (t e ) = —m H N 0 —- i r c [t e ) 

o 1 + OL 


(C3) 


Equations (Cl) and (C3) can be solved for t e , and 


te. = t, 1 


1 

rjco 


3^ + r} c0 - 1 


Zu + 2(2u 2 + or) /(I + a) 


1/2 


(C4) 


At times t > t e , clump expansion is halted, and the clump is now confined by the pressure at the 
base of the PDR flow. The clump slowly loses mass, and shrinks to eventually become completely 
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photoevaporated. The time evolution of clump mass and size can be obtained by differentiating 
equation (C3) which also holds for t > t„ with respect to time and equating this with the mass 

rate (equation Bll ). Thus, 


1 + a 


;(l > te) = r e (te) - 2 CpDR 2 V 2 + oS l te ^ 


2(2;/ 2 + a) / r c (t) 


m c (t > t.) = m c0 ^ (1 + a) y rd) 

The photoevaporation timescale ,s obtainetl by setting the clump radius to zero at t rE to give 


tpE ~ 


2v 2 + a\ r c (0) + c c t e 


1 + a 


3cpdr 


+ tet 


loss 

(C5) 

(C6) 

(C7) 


and t e is given by equation (C4). 


D. Numerical hydrodynamics code 

The fluid equations of motion for the system are solved using a 1-D spherical Lagrangian 
hydrodynamics code. The equations are solved using a finite difference method with a numerical 
viscosity term, added as a pseudopressure, for accurate handling of shocks and discontinu tie 
(Richtmyer k Morton 1967, Bowers k Wilson 1991). An isothermal equation of state 1S ' lsed 
throughout. The accuracy and stability of the scheme was checked with standard test problems with 
known solutions. The computational grid has 1000 equally spaced radial zones, except for the central 
20%, which was spaced logarithmically in radius. This was done to increase spatial resolution at the 
centre, and it provides a more accurate calculation of any shock-induced collapse. 

The FfJV field is gradually turned on over a timescale, t. F uv, which is varied to accommodate both 
the impulsive and slow heating cases described in §4. The temperature of an outer shell of material is 
thus raised continuously from T, to T FD r, reaching its maximum value of the PDR temperature at a 
time truv The temperature T exponentially drops off with the square of the column, N, into the 
cloud TIN) = T + (TpnR - T c fe-f"/""! . This profile was chosen to closely mimic the T dependence 
with column predicted by PDR models for the values of Go and n under consideration (Kaufman et al. 

1999). 



Table 1: List of symbols frequently used in this paper 


Symbol Meaning 

c c Thermal sound speed in cold clump gas 

cpdr Sound speed of heated, warm PDR gas 
Go Strength of the FUV field in Habing units, 1.6 x 10 -3 erg s _1 cm -2 
No Fiducial column density equal to 10 -21 cm -2 
rico Initial gas number density of clump 
m c o Initial mass of clump 

r c0 Initial radius of clump 

t c Sound crossing time in cold clump gas, r c0 /c c 
t e Expansion timescale for clumps with 7j c0 > rjcru 
tFuv Timescale on which the cold clump gas gets heated to the maximum temperature 
in a given FUV field 
tpE Photoevaporation timescale 

t s Shock-compression timescale for clumps with r ] c0 < r] crit 
Vf Velocity of the Ionization Front moving into the GMC 
Npdr Thickness of the PDR in the interclump medium of the GMC, bounded by the IF and DF 
a Initial ratio of turbulent to thermal pressures in the cold clump gas 

(3 Initial ratio of magnetic to thermal pressures in the cold clump gas 

7 Power law exponent of density varying with magnetic pressure 

5 0 Initial thickness of PDR shell on the surface of a clump (= No/n c o) 

Tj cr Critical column density for complete photoevaporation of pressure-confined clumps 
r ko Normalized initial column density to centre of clump, n c 0 r c o/No 

Vcrit Critical column density in a given FUV field towards which unconfined clumps evolve 

rjj Final column density of pressure-confined clumps in an FUV field 
A Photoevaporation parameter; A = 1 for 77 = r) crit 
u Ratio of sound speeds in warm and cold gas, cpppjc c 
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Fig. 1.— A schematic diagram of PDRs on large scales (e.g. Rosette) and small scales (e.g. Orion 
Bar). The figure shows blister H II regions formed by the massive stars, and the PDR bounded by the 
ionization and dissociation fronts. Clumps in the PDR are exposed to FUV radiation and lose mass by 
photodissociation and/or photoevaporation. 
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L °g v c0 

Fig. 2. — Logarithmic parameter plot of the ratio 77 of the column density through the clump to the 
fiducial column N 0 , and the ratio u of the sound speed in the heated PDR gas to that in the cold clump 
gas, for impulsively heated clumps. The critical column density r) cr u is plotted against u for different 
values of a and /3, with 7 = 4/3. Also shown are representative trajectories of evolution for clumps with 
initial column densities less than, and greater than, T] cr i t . A clump with r] c o < r] cr it is compressed by 
a shock to rj crit and moves right on this parameter plot and then proceeds to photoevaporate steadily. 
A clump with i] c o > rj CTlt expands so that its column density becomes equal to ri CTlt and continues to 
photoevaporate. The initial turn-on of the FUV field (dotted line) moves the column slightly to the left 
because a column jV 0 is quickly dissociated and lost from the clump. Gravity is ignored in this figure 
(see Figure 5). 




- 35 - 



Fig. 3. — Schematic diagram of evolution of shock-compressed clump. The cold clump is indicated by 
the lightly shaded region, the dotted region represents the heated warm PDR surface of the clump, and 
the dark shaded region shows shock-compressed gas. Panel c shows the propagation of the shock front 
into the clump and panel d depicts the evolution of the clump after being compressed by the shock into 
a small dense core. 
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Fig. 4. — Time evolution of radius (panel a) and mass (panel b) of cold clump gas for a clump with 
rjco = 10, v = 5 that undergoes shock compression. The mass and radius are given as ratios to the initial 
mass’ and radius, and time is in units of sound crossing times in the initial cold clump. The solid line 
is the analytical result, showing an initial rapid decrease in clump radius followed by a slower phase of 
clump radius evolution. The dashed line is the result from the 1-D hydrocode and is seen to agree well 
with the analytical curve. 
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p ig 5 — This figure shows the rj - v parameter space and demarcates the regions where gravitational 
collapse can be triggered by shock compression (for 7 = 1). The dashed line indicates the co umn 
densities of clumpsJwith thermal pressure 10 6 cm" 3 K, which are initially unstable to collapse Clumps 
with this pressure and columns to the left of the line are initially stable. The regions to the right of the 
solid line are where clumps are driven to collapse by shocks induced by photoevaporation, for initial 
thermal pressures in the clump P T = 10 6 . For higher pressures P c0 , both the dotted (r? a P c0 ) and 

solid lines (7/ oc P c 2 0 /3 ) move to the right. 



- 38 - 




t/t c 


Fig. 6 — Time evolution of radius and mass for an initially expanding clump, with t? c0 = 100, u = 3. 
Radius and mass are given as ratios of their initial values and time is in sound crossing time units, t c . 
The figure shows the analytical result (solid line) and that from the numerical code (dashed line). The 
clump radius increases in the beginning as the clump expands and then shrinks as the clump loses mass 
due to photoevaporation. 
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Log 77, 

Fig. 7. — The evolution of column density r]j of the clump against the parameter v. As soon as the clump 
is exposed to an FUV field, an outer column Nq is instantaneously heated, and thus for u marginally 
greater than 1, at t = 0, the final column density 7// = r) c0 - 1. The figure shows how the final column 
density of the equilibrium clump configuration decreases 'with an increasing FUV field. Also shown are 
the contours depicting the percentage of mass lost from the clump as it photoevaporates. 
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Fig. 8. — Density, velocity, pressure and temperature as a function of radius at times t = 0.05,0.1, 
and 0.1 5t c . The density, pressure and temperature are scaled to the initial values in the clump and the 
velocity is in units of the clump thermal speed c c . 
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Fig. 9. — This figure shows the mass in the clump, core and warm PDR shell as a function of time. 
At collapse, the clump mass drops sharply to zero. Also shown in dotted lines are the clump and PDR 
shell masses for the case where 7 = 4/3. 



Fig. 10 — This figure shows the evolution of mass with time for a slowly heated, pressure-confined 
clump (Bonnor-Ebert sphere) with varying central to surface density ratios, x = 1,3 and 6. In this case, 
T) c0 = 10, v = 5, and t FUV = 5 t c . The corresponding analytical solutions for a constant density sphere 
(dotted line) and an isothermal sphere (dashed line) are also shown for comparison. 
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Fie 11 - Observational data for some PDRs on the rj - u parameter plot for the turbulent (panel a) 
and pressure-confined (panel b) clump models. There is a slight indication that the columns he along 
the lines expected in the turbulent model of clumps, with the exception of NGC 2023. 
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